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ABSTRACT 

We present the results of a thorough study of the nucleus of the Sgr dwarf spheroidal galaxy (Sgr 
dSph) and of the bright globular cluster M54 (NGC 6715) that resides within the same nucleus 
(Sgr,N). We have obtained accurate radial velocities and metallicity estimates for 1152 candidate 
Red Giant Branch stars of Sgr and M54 lying within ^ 9' from the center of the galaxy, from 
Kcck/DEIMOS and VLT/FLAMES spectra of the infrared Calcium II triplet. Using both velocity 
and metallicity information we selected two samples of 425 and 321 very-likely members of M54 and 
of Sgr,N, respectively. The two considered systems display significantly different velocity dispersion 
profiles: M54 has a steeply decreasing profile from r = 0', where a ~ 14.2 km/s, to r ~ 3f5 where 
it reaches a ~ 5.3 km/s, then it appears to rise again to cr ~ 10 km/s at r '--^ 7'. In contrast Sgr,N 
has a uniformly flat profile at cr ~ 9.6 km/s over the whole 0' < r < 9' range. Using data from the 
literature we show that the velocity dispersion of Sgr remains constant at least out to r ~ 100' and 
there is no sign of the transition between the outer flat-luminosity-profile core and the inner nucleus 
in the velocity profile. These results - together with a re-analysis of the Surface Brightness profile of 
Sgr,N and a suite of dedicated N-body simulations - provide very strong support for the hypothesis 
that the nucleus of Sgr formed independently of M54, which probably plunged to its present position, 
coincident with Sgr,N, because of significant decay of the original orbit due to dynamical friction. 
Subject headings: galaxies: dwarf — globular clusters: individual(NGC 6715) — stars: kinematics — 
galaxies: nuclei — galaxies: individual (Sgr dSph) 
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1. INTRODUCTION 

The occurrence of stellar nuclei at the photomet- 
ric center of several dwarf elliptical galaxies has been 
the subject of s tatistical investigations since the pio - 
neer ing study bvlBinggeli. Sandage fc TainmannI ()1985D 
and iBinggeli. Tammann fc Sandagd (|1987f ). and nucle- 
ated dwarf ellipticals (dE,N) have become a generally 
recognized and wel l studied c lass of galaxies in their own 
right (see .Gallagher fc Wvsa ,1994: Ferguson fc Binggehl 
11994 and references therein). In the last few years the 
results from systematic studies performed with the in- 
strumentation on board of the Hubble Space Telescope 
(HST), has revolutionized the field, providing strong ev- 
idence supporting the possibility of an intimate connec- 
tion between the process of nucleation and the process 
of galaxy formation as a whole. In particular: 

• It was generally accepted that stellar nuclei 
occurred in a significant but minor frac- 
tion of dwarf elliptical g a laxies (~ 25% 
IBinggeli. Sandage fc TammannI |1985[ ): on the 
other hand, the high resolution HST analysis 
of a large and well selected sample of Virgo 
dEs by ,C6te et a l. (2003, hereafter C06) has 
revealed that the fraction of dE,Ns (/„) can be 
as high as 66% < fn < 87%. Strong support 
to a significant upward revision of /„ comes 
also from the recent gro und-based study by 
iGrant. Kuipers. fc Phillippi (j2005s). Hence, nucle- 
ation appears as the natural status of dEs instead 
of an exceptional occurrence (see C06 for details 
and discussion). 



An HST study o f a lar^ e sample of Sa-Sd galaxies 
by iBoker et a l. (2004) found that a similarly 
high fraction of stellar nuclei is found also among 
these late type spirals, /„ ~ 77%. In general, the 
most recent studies agree in finding a fraction of 
nucleation larger than 5 0% in any kind of galaxy 
(see, for example |Carollo. Stiavelli fc Maclj 
199S: [Balcells. Graham fc Peletieil l2007l : 
Ferrarese fc Cotel l2006t iBoked l200'7r and ref- 



formation, while its further evolution into a SBH or 
a stellar nucleus would depend on various "local" 
circumstances. 



erences therein), with the only exception of those 
brigh ter than Mb — —20.3 mag ()Ferrarese fc Cotel 
l2006l hereafter FC06). Hence, nucleation seems 
to occur very fre quently in both dwarf an d gi- 
ant galaxies (see iGraham fc Guzmiiil l2003l and 
references therein, for a deeper discussion of the 
dwarf/giant dichotomy). 

• Stellar galactic nuclei are found to obey the 
same scaling relation that links Supermassive 
Black Holes (SBH) w ith their host galaxies 
(iWehner fc Harris! [200l C06, FC06). It has been 
suggested that ^^...SBHs and stellar nuclei are 
nothing but complementary incarnation of Central 
Massive Objects - they likely share a common for- 
mation mechanism and follow a similar evolution- 
ary path..." (FC06). In this view, the accumula- 
tion of a compact overdensity of baryons at the very 
center would be a common process during galaxy 
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• The similarity between stellar nuclei, the brightest 
globular clusters and the recentl y discovered Ultra 
Comp act Dwarf (UCD) galaxies (jPrinkwater et ahl 
l2003f ) clearly suggests possible links among 
these classes of s tellar systems (jFreemani 19931 : 



Lotz et al.' '2001; 'Lotz, Miller fc Ferguson' '200 
Bass i no, Muz zio & RaboUi 1994; Hasegan et al 
20051 : iFederici et all 12007) : iBokCTi l2007l) . " In par 



ticular it has been suggested either that nuclei 
can f orm by merging of pre-existing globular clus- 
ters (jLotz et al.ll200lL and references therein), or 
that some objects currently classified as globular 
clusters may be in fact galactic nuclei, the only 
relics of the tidal shredding of their host galax- 
ies t hat once orbited the Milky Way and/or M31 
(seelMack ev fc van den Bergh''2 005l : IFederici et al.l 
120071 : iBrodie & Strader 2006, and references 
therein), or even that most/all globulars were nu- 
clei of Galactic build i ng blocks (i.e. dwarf galaxies, 
see Freema,nl Il993t iFreeman fc Bland-Hawthorn! 
I2OO2I : lB6keij|2007l . BS. and references therein). 



Within the scenario briefly outlined above, the cru- 
cial relevance of any independent observational fact that 
may provide some insight on the process of nucleation is 
clearly manifest. The mechanisms for the formation of 
stellar nuclei that have been suggested until now belong 
to two broad classes: (a) the orbital decay and merg- 
ing at the center of the parent galaxies of pre-existing 
stellar systems (star clusters), and (b) the accumula- 
tion of gas (of various origins) in the very center of a 
galaxy and its subsequent transfo rmation into stars (see 
iGrant. Kuipers. fc P hillipps 2005, for a detailed descrip- 
tion of various proposed flavors of these classes of mod- 
els). 

At present, all the observational constraints we have 
in hand come from the study of the integrated proper- 
ties of samples of distant galaxies/nuclei. The nearest 
galaxies that are known to host a stellar nucleus are 
M32, M33 and NGC205, in the Local Group. Even in 
the most favorable case of NGC205, the structure of 
the nucleus can be studied in some detail but it can 
be resolved into indiv idual stars onl y partially, even 
with HST imaging (Valluri et al.!l2005f ). In this context 
the c ase of the Sagittarius dwarf sp heroidal galaxy (Sgr 
dSph, llbata. Irwin fc Gilmorelll994f) constitutes an abso- 
lute unicum ( Monaco et al Jl2005a|, hereafter M05a). 

The Sgr dSph is a dwarf sa tellite (L - 2 - 5 x IO'^Lq) 
of the Milky Way (jibata et a l. 1997) whose main body 
is at Dq — 26.3 ± 1.8 kpc fro m us and Rgg — 18 .7 kpc 
from the cente r of the Galaxv ([Monaco "etal][200l . It IS 
devoid of gas (jBurton fc Lockmanl 11999! ) and its stellar 
content is dominated by an intermediate-age (~ 6— 8 Gyr 
old, see jBellazzini et al. 200_6a„ and references therein) 
metal-rich ((Z) ~ § Z0. lMonaco et al.ll2005"bl and refer- 
ences therein) population, w hile :^ 10% of old and metal 
poor stars are also present ([Monaco "eral][200l . The 
Sgr dSph is currently beijng destroyed by the Galactic 
tidal field (jibata et al.![l997l) . The tidally stripped stars 
form a huge and coherent tidal tail system that can be 



The nucleus of the Sagittarius dwarf galaxy and M54 
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obs erved over the w h ole celestial sphere (S gr Stream, 
see llbata et alJ 120011 : iMaiewski etaH l2003l and refer- 
ences therein). There are four globular clust ers within 
the main body of the galaxy (llbata, Irwin fc Gilmor^ 
1994 iDa Costa fc Armandroi 119951 : iMontegriffo et aTl 
19981 ) and probably several others that have been 
stripped from the m ain body and are now associated 
with the Sgr Stre a m dBellazzini. Ferraro fc Iba ta"2003al: 
Bellazzini et al. ] I2003bt IC arraro. Zinn fc M oni Bidinj 
20071) . 

Particularly n oticeable, how ever, is the very bright 
(My = -10.0. iHarrisI Il99^ metal-poor (\F e/H\ = 
-1.55) (jBrown. Wallerstein fc Gonzalej Il999l ) globular 
cluster M54 (NGC6715), located exactly at the center of 
the light distribution of the Sgr galaxy. For this rea- 
son, immediately after the discovery of the galaxy it 
was sugge sted that M54 could be the actual nucleus of 
Sgr d Sph (H assino fc Muzzioiri995l : ISaraiedini fc LavdenI 
119951). Some doubt on thi s hypo thesis was advanced 
by ' Da Costa fc Armandrofil (|1995f ) based on the cir- 
cumstance that the integrated color of M54 is much 
bluer than that of its host galaxy, while the oppo- 
site was generally be l ieved to occur in dE,Ns; however 
ISaraiedini fc LavdenI ()1997l ) noticed that there were ex- 
ceptions to this rule. In fact, the latest studies suggest 
that the opposite may be true (Lotz. Miller & Ferguson 
I200I see also C06). Since the mean metallicity of M54 
and of the population dominating the Sgr galaxy differ by 
one full dex, their Red Giant Branch (RGB) stars (and 
in some cases also their Horizontal Branch stars, HB) 
can be easily discriminated from each other in a Color 
Magnitude Diagram (CMD); at the same magnitude the 
RGB of Sgr is much redder than that of M54. Se l ecting 
gen uine Sgr stars in t his wa v lLavden fc Sarajedinil (j2000l ) 
and lMaiewski et"al] (|2003D independently found an over- 
density of Sgr stars that appear concentric with M54 
and have a similar spatial scale. Using the same tech- 
nique and a very large optical photometry database M05a 
demonstrated that Sgr has actually a nucleus of metal 
rich stars (in the following we will refer to this structure 
as the nucleus of the Sagittarius galaxy, Sgr, N for brevity, 
if not otherwise stated). This stellar structure has the 
same center as M54, to within the uncertainties, but it 
displays a different surface brightness profile, suggesting 
a different origin from the cluster. The strong incom- 
pleteness in the innermost ~ 10" prevented M05a from 
obtaining accurate estimates of the structural parameters 
of the Sgr nucleus: they derived —10.0 < My < —7.6, 
M v(0) 5 . 18.5 mag/arcsec^, and, with a tentative fit of 
a lKingr()1966l ) model < 0.21' and rt ~ 17' (where 
C=log(r(/rc)~ 1.9; we will refine these parameters in 
Sect. 2, below). They also found that the observed prop- 
erties of the nucleus of the Sgr galaxy were fully compati- 
ble with those of known n uclei of dwar f ellipt icals. More- 
over, the recent study bv lSiegel et al.l ()2007| ) has shown 
that while the stellar budget of Sgr,N is dominated by 
the same intermediate- age population found in the sur- 
rounding galaxy, there is clear evidence of more recent 
and repeated episod es of star formation (no t observed in 
extra- nuclear fields. [Bellazzini et al.|[2006af ). pointing to 
a history of subsequent phases of re-accumulation of gas 
followed by sta r formation bursts , very similar to what 
is observed by iRossa et al.l ()2006D in a large sample of 
distant galactic nuclei. 



In any case, the key result of M05a is that Sgr is a 
nucleated galaxy independently of the presence of M54; 
even if one were able to remove the cluster from the 
galaxy by magic, a nucleus made of typical Sgr stars 
would still be there. To explain the strict spatial coinci- 
dence of M54 and the nucleus, M05a proposed two possi- 
bilities: (a) the nucleus formed in situ from processed Sgr 
gas that has fallen to the bottom of the galactic poten- 
tial well, and M54 was (independently) driven into the 
same place by dynamical friction; (b) M54 was born in 
the present place (or was driven there as in case a) and 
formed a co-spatial overdensity of Sgr stars by capturing 
them (or the gas from which they formed) from the sur- 
roundings. The key difference between the two hypothe- 
sis can be summarized by the questions "did M54 provide 
the mass seed around which the metal-rich nucleus of Sgr 
was assembled? is M54 the main contributor of the mass 
budget in the innermost ~ 100 pc of the Sgr galaxy?", 
that enclose the basic point of hypothesis (b) . 

Hypothesis (a) was found to be fully compatible with 
all the available information, while it was not possible 
to consider case (b) i n more depth. In a recent paper, 
iMieske fc Baumgardtl ()2007l ) studied in detail the pro- 
cess of the capture of "field" stars from an intervening 
star cluster and they concluded that, even in the most 
favourable cases, the fraction of stars captured in a Hub- 
ble time is negligible, less that 10 ~^ of the cluster mass. 
The ra nge of conditions studied bv lMieske fc Baumgardtl 
(|2007t l includes the "M54 within Sgr" case, hence their 
conclusions are fully applicable here. Nevertheless, the 
results by Micskc fc Baumgardt (2007_) do not exclude 
the possibility that M54 operates as a collector of pro- 
cessed gas at the center of Sgr. 

However it seems quite reasonable to expect that if 
case (b) is true, the stars of the Sgr,N structure should 
share the same kinematics as M54 stars, or, at least, their 
kinematics should be compatible with a mass-foUows- 
light model, or, finally, their kinematics should be differ- 
ent from that of extra-nuclear Sgr stars, as they would 
be bound to M54 a nd would orbit within its potential 
(iGilmore et al.ll2007f ). Note that M54 and Sgr,N are en- 
closed within the innermost 10' of a galaxy whose core ra- 
dius and limiting radius are as larg e as rc = 224'zfclO' an d 
rt = 1801' ± 112', respectively (M aiewski e^a^l2003^ . 
Hence, even if the main body of Sgr is clearly undergo- 
ing tidal disruption the kinematics of the nuclear region 
should be unaffected and is consequently expected to pro- 
yide a trustworthy insight on the inne r mass distribution 
(|Munoz. Maiewski fc Johnston|[2007f) . 

To follow-up these ideas and to obtain a deeper insight 
into this unique case of a galactic nucleus that can be 
studied in vivo on a star by star basis, we performed 
an extensive spectroscopic survey, principally aimed at 
the study of the kinematics of Sgr,N and M54. Here we 
present the main results of this survey. The plan of the 
paper is as follows. As a preliminary step we re-analyse 
the Surface Brightness profile of Sgr,N with new data 
in Sect. 2. In Sect. 3 we describe the observation and 
data reduction of our spectroscopic survey, we test the 
reliability and the accuracy of our measures, the observed 
metallicity distribution is briefly discussed and, finally, 
the selection of the samples is described in detail. Sect. 4 
is devoted to the kinematical analysis of the M54 and 
Sgr,N samples and in Sect. 5 we briefly present the results 
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F808W-FB14W F808W-FfiI4W 

Fig. 1. — Upper left panel: the annuli and semi-annuli adopted 
for radial star counts are plotted on the X,Y map of the stars 
having F814W< 20.0; X,Y are the local coordinates as defined in 
Sect. 3., below, but in arcseconds. Upper right panel: overall 
CMD; note the difference among the steep and star rich RGB of 
M54 and the redder and sparser RGB of Sgr.N, for F814W< 17.0. 
Lower panels: the selection of M54 and Sgr,N RGB stars adopted as 
density tracers is illustrated for two radial ranges. The horizontal 
lines enclose the magnitude range of the selection, the bluest and 
the reddest red lines select the RGB stars of both systems and 
the middle red line separates the sample into M54 members (bluer 
than the line) and Sgr,N members (redder than the line). 

of a suite of N-body simulations aimed at the study of 
the orbital decay of M54 within the Sgr galaxy. Finally, 
in Sect. 6 we summarize and discuss our results. 



2. REFINING THE STRUCTURAL PARAMETERS OF 
SGR.N 

To gain a better insight into the Surface Brightness 
(SB) profile of Sgr,N, we reduced archive data obtained 
with the Wide Field Channel (WFC) of the Advanced 
Camera for Survey (ACS) on board of the Hubble Space 
Telescope (HST). This is a set of short and long ex- 
posure F814W and F606W images of the center of 
M54, taken within th e GO 10775 Treasury Program 
(|Saraiedini et al.|[2007f) . The relative and absolute pho- 
tometry of individual stars has been obtained following 
the same steps described in lSollima et al.l (2007) for data 
taken from the same survey. Further details and the over- 
all analysis of these data will be presented elsewhere; an 
independent study of the stellar populations of M54 and 
Sgr.N from the same images has been recently presented 
bv ISiegel et al.l ([20071) ■ In the present context, we limit 
ourselves to complementing the SB profile of Sgr,N ob- 
tained by M05a with a couple of points in the innermost 
10", using the same methodology as M05a, and, as a con- 
sequence, to obtain stronger constraints on the structural 
parameters of Sgr,N. 

Fig. [1] clearly shows that the new ACS photometry al- 
lows one to discriminate the RGBs of M54 and Sgr,N by 
color even in the innermost 10". Using the selection il- 
lustrated in Fig. [T] we were able to reliably measure the 
ratio of the density of tracer stars (RGB and AGB, in 
the present case) in the two samples in any given ra- 
dial annulus psgr,N / Pmsa, as it scales as the ratio of the 



TABLE 1 

Surface Brightness profile of the innermost 60" of Sgr 
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Note. — and r j are the limits of the bins, Vm is the average 
radius of the sample. 



TABLE 2 

Parameters of the King model that best fits the Sgr,N 

PROFILE. 



Parameter Estimate Uncertainty Units 





0.05 


± 0.01 


arcmin 


c 


1.90 






rt 


4.00 


± 0.8 


arcmin 




0.42 


± 0.08 


arcmin 


ri 


~10.5 




arcmin 




15.3 


± 0.2 


mag/arcsec-^ 


My 


-7.8 


± 0.2 


mag 



Note. — rt is the tidal radius of the same King model, while 
r; is the observed limiting radius, defined as the radius at which 
the observed profile of the nucleus appears to join the profile of the 
overall Sgr galaxy. All the reported values and uncertainties have 
been estimated by assuming the reported value of C. The reported 
values of the surface brigthness are not corrected for inters tellar 
extinction. is the half-mass radius of the best fitting K ind 

(1966) model: it can be considered equivalent to the half-light 
radius, in the present case. The half-light radius of the lKingI 11196^) 
best fitting the SB profile of Sgr,N is = 0.23'. 

number of stars per annulus Nsgr,N /NMbA- To avoid 
the inter-chip gap we used semi-annuli for R > 10". 
The quantity 2.b\og[Nsgr,N /NMbi)^ in turn, scales as 
the difference in Surface Brightnes s between the two 
syste n is (see M05a for details, a nd iRenzini fc Buzzonil 
119861 : IRenzini fc Fusi Peccil Il988l . for the basis of the 
underlying evolutionary flux theorem)'^. No attempt 
was made to estimate the ellipticity of the system; as 
the system appears quite spherical at large radii (see 
M05a) we assumed spherical symmetry. The coordinates 
of the centers of M54 and Sgr,N, as derived from the 
2(7 clipped average X and Y of stars selected by color 
as in Fig. [H are found to coincide to within < 2", 
in excellent agreement with M05a. The effects of the 
radial variation of completeness on the density ratio 
should be practically null, as we selected the M54 and 
Sgr,N samples in the same (bright) magnitude range 
(13.8 < F%UW < 16.4). To convert the estimated SB 
differences into an absolute scale we adopted the best- 
fit mo del profile of M54 by iTrager. King fc Diorgovskil 
(|1995f ) as a "zero-point" , consistent with M05a. In prac- 
tice, for each observed Aii^3^^^'^^^{Ri) we obtained the 
surface brightness of Sgr,N at radius Ri, fiy^^{Ri), as 
^^f''(i?,) = p^^^^iR,) + Afi^a-'-^^^iR,). The derived 

We used s ynthetic CMDs produced with the dedicated BASTI 
web interface (|Cordie r et al.| [ 2007D to check that a population of 
age = 12 Gyr and [Fe/H]=— 1.5 (taken as a proxy for M54) and a 
population of age = 6 Gyr and [Fe/H]=— 0.4 (taken as a proxy for 
Sgr,N) having the same total V luminosity, places (approximately) 
the same number of stars in the selection windows shown in Fig.[T] 
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O M54 (MOSa) 
A Sgr (M05a) 
▲ Sgr (ACS) 



1 2 3 

Log r(arsec) 



Fig. 2. — Surface brightness profiles of M54 (pentagons) and of 
Sgr (triangles). Continuous open symbols with errorbars are from 
MOSa. The filled triangles have been obtained from the density 
ratio computed here from HST/ACS data and are normalized by 
adding them to the corresponding SB of the M54 best-fit model, 
here marked by dotted pentagons. The long dashed line is the best- 
fit King (1966') model found by Trager. King & Djorgovski (1995) 
for M54; the sho rt dashed line is the best-fit King model found by 
IMaiewski et al.l (|2003) for the overall main body of the Sgr galaxy; 
the dotted line is the adopted best-fit model for Sgr,N and the 
continuous thick line is the sum of the last two models. 

portion of the profile is reported in Tab. [T] and plotted in 
Fig. [21 together with the portions of the M54 and Sgr,N 
profiles derived by MOSa. 

The SB of Sgr,N is estimated at five radii. The three 
points at i? > 10" nicely overlap with the observations 
by MOSa, showing that the two profiles are in good agree- 
ment. The two points at i? < 10" provide very strong 
new constraints on the fit of the overall profile. A C=1.90 
King model with /iv'(O) = IS. 4 and = 3.0" provides 
a good fit of the profile from the innermost observed 
point, at i? ~ 3", to i? - 120" = 2'; for i? > 2' the 
observed profile shows a much gentler decline with re- 
spect to the model, joining the fiat profile of the core of 



Sgr a.t R — ri ~ 10.5' instead of R 



4.0'. While 



there is no particular physical reason in adopting a King 
model to represent Sgr,N, it provides an easy and satis- 
fying way to parametrize its inner SB profile. Moreover, 
while the best-fit King model fails to reproduce the ob- 
served profile in the range 2' < i? < 10. S', the fraction of 
the total Sgr,N light enclosed in this range amounts to 
just a few per cent; therefore the adopted model provides 
a reasonable description of the distribution of the bulk 
of the system light. 

The derived parameters are listed in Tab. O Our re- 
sults are fully compatible with the limits and the edu- 
cated guesses by MOSa; the most remarkable difference 
is in Tc that was tentatively guessed as larger than that 
of MS4 by MOSa and, in fact, is found to be signifi- 
cantly smaller here. The absolute integrated V mag- 
nitude (My — —7.8) is just slightly brighter than the 



upper limit by MOS (My 



< 



-7.6). Assuming the best-fit 



models as correct, Sgr,N has a central V luminosity sur- 
face density that is slightly more than half that of MS4, 
its total V luminosity is j of the cluster and, finally both 
its core and half-light radii are ~ one half of those of the 



cluster"'. 

With these new and much tighter constraints on the 
light distribution of Sgr,N in hand, we will proceed in 
the following sections to perform a comparative study of 
its kinematics and those of MS4. 



3. THE SPECTROSCOPIC SURVEY: OBSERVATIONS, DATA 
REDUCTION AND SELECTION OF THE SAMPLES 

The spectroscopic observations on which this study 
is based were obtained with the multi-object spectro- 
graphs DEIMOS on the Keck 2 telescope and FLAMES 
on the VLT-UT2 telescope, as listed in Tab. El The 
DEIMOS observations were undertaken using two differ- 
ent observing strategies: the first used the standard "slit- 
let" mode (the slitmask app roach used by the DEEP2 
team, e.g.. lDavis et al.ll2003f ). where short slits of mini- 
mum length 4" and width 0".7 were milled; the second 
mode used smaller long and 0'^7 wide "holes" to 
allow larger multiplexing. In both cases, the high res- 
olution 1200 lines/mm grating was employed with the 
OGSSO blocking filter, giving a resolution of approxi- 
mately 1.4l FWHM. The slitlet spectra covered 6500- 
9000A and were extracted and wavelength calibrated us- 
ing the DEEP2 pipeline software; in contrast, the "holes" 
spectra wer e extracted using software developed by our 
own group (jlbata et al.ll200Sl ). with the final extracted 
spectra covering only the region 8400-87S0A (purely due 
to a limitation in the software). 

For the DEIMOS observations we selected candidate 
MS4/Sgr,N RGB and Red Clump (RC) stars lying within 
i? < 9' from th e center of the s ystera s from the wide-field 
photometry of [Monaco et al.l (|2002f ). A limiting magni- 
tude of I = 18 was adopted, which ensures that we probe 
the Red Giant Branch to approximately one magnitude 
below the Red Clump (see Fig. [3]) . 

Additional FLAMES observations were downloaded 
from the ESO archive; these targeted almost exclusively 
MS4 stars, with only a small contamination from other 
populations. The FLAMES spectra were obtained with 
the high resolution setting HR21, which covers the cal- 
cium triplet from 8484-87S7A, with a resolution of O.SA 
FWHM. The FLAMES data were extracted and cali- 
brated using the "girBLDRS" pipeline^, developed for 
the European Southern Observatory by Geneva Obser- 
vatory. Being of higher resolution, the FLAMES observa- 
tions can be used to assess the accuracy of the DEIMOS 
data, but more importantly FLAMES is a fiber-fed spec- 
trograph and is therefore immune to the velocity errors 
that can arise from mis-centering of slit spectrographs 



* r^j is the half-mass radius of the best fitting lKintj (119661 ) model: 
it can be considered equivalent to the half-light radius, in the 
present case. Note that using the same definition, the half-light 
radius of M54 is = .79', significantly larger than the value re- 
ported bv IHarrTi I I1996I . = 0.49'). The latter number is very 
similar to what we obtain for the rjj of the KinJ I I1962I) empirical 
model that best-fits the profile of M54 (having the same C and 
Tc adopted here), i.e. = 0.47'. For C> 1. dynamical IKin j 
(jT966) models have rj, larger than IKintj 119621) empirical models 
by more than 10%; the difference becomes larg er than a factor of 
two for C > 2.0. The half-light radius of the iKmg (1962) best 
fitting the SB profile of Sgr,N is = 0.23'. Hence independently 
of the adopted family of models, the half-light radius of M54 is a 
factor of 2 larger than Sgr,N one. Preliminary tests suggest that 
available compilations of rj, likely include quantities estimated in 
various non-homogeneous ways. 

^ See jhttp : //girbldr s ■ sour cef orge ■ netT] 
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TABLE 3 
Spectroscopic observations. 



Field 



Instrument 



Date Exposure 



1 


FLAMES Medusa 1 28/07/2005 


3 


X 


845 s 


2 


FLAMES Medusa 1 28/07/2005 


3 


X 


845 s 


3 


FLAMES Medusa 2 28/07/2005 


3 


X 


845 s 


A 

4t 


FLAMES Medusa 1 


20/08/2005 


q 
o 


X 


/yo s 


5 


DEIMOS holes 


31/08/2005 


3 


X 


900 s 


6 


DEIMOS holes 


31/08/2005 


3 


X 


900 s 


7 


DEIMOS holes 


02/10/2005 


3 


X 


950 s 


8 


DEIMOS holes 


02/10/2005 


3 


X 


950 s 


9 


DEIMOS holes 


03/10/2005 


3 


X 


720 s 


10 


DEIMOS holes 


27/05/2006 


3 


X 


300 s 


11 


DEIMOS holes 


27/05/2006 


3 


X 


300 s 


12 


DEIMOS slitlets 


28/05/2006 


3 


X 


720 s 


13 


DEIMOS slitlets 


23/09/2006 


3 


X 


300 s 


14 


DEIMOS sUtlets 


25/09/2006 


3 


X 


300 s 



such as DEIMOS. 

The radial velocities of the target stars were obtained 
by cross correlating the observed spectra against an ar- 
tificial template of the Ca II lines, using th e same ap- 
proach described and discussed in detail in llbata et al.l 
(|2005D and iBattagUa et al. I (|2008f ). The determination 
of the uncertainty on the velocity is based on a weighted 
sum of Gaussian fit errors to the individual CaT lines 
as described in detail in the thorough discussion by 
iBattaglia et al. I ()2008i) . who also demonstrated that the 
velocity errors computed in this way are reliable. An 
additional sanity check was provided by the comparison 
with the scatter of the velocity measurements derived 
from the thre e individual C a II lines separately (see, 
again, Battagha et al. II2008D . 

The metallicity of the target stars were estimated from 
the combined equivalent width of the CaT lines. To 
this end we i r nplern ented the calibration described by 
iCarrera et al.l ()2007[ ). which it is claimed works well up 
to high [Fe/H] values. As usual for CaT calibrations, 
[Fe/H] is function of the combined equivalent width and 
of the difference between the V magnitude of the star 
and the mean level of the HB /R R Lyrae of the system 

V — Vhb- For full consistency with'C arrera et ahl (j2007[ ). 
the equivalent widths of the CaT lines were measured 
according to their prescript ions. The metallicit y mea, - 
surements were put on the ICarretta fc Grattonl ()1997D 
scale, ass uming Vhb = (Vr.b.Lv) = 18 .17ib0.01 and E(B- 
V)=0.14 (jLavden fc Saraiedinil[2?)0nh . 

3.1. Photometric properties of the sample 

From these data we constructed a final catalogue con- 
taining the radial velocity and the metallicity of 1152 
stars in the innermost 15' x 10' of the Sgr galaxy, as 
displayed in the lower panel of Fig. [31 The CMD of all 
these stars is shown in the upper panel of the same fig- 
ure. The RGB sequences of the two systems are quite 
well separated in the CMD down to the faintest stars in- 
cluded in our sample (/ ~ 18). However, we also selected 
stars lying on the Red Clump of the Sgr,N population, at 

V — I ~ 1.15 and 16.8 < I < 17.3: this feature overlaps 
with the RGB of M54, hence we will need some further 
information to establish the membership of these stars, 
in addition to their color, magnitude and radial veloc- 
ity. In particular, the metallicity would be very useful 
to disentangle the two populations in this range. M54 



has \ Fe/H] — —1.55 (iBrown. Wallerstein fc Gonzaled 
119991 : iDa Costa fc Armandroail995D. w ith a small dis- 
persion fsee ISaraiedini fc LavdenI Il995l and Sect. 3.4, 
below). In contrast, the metallicity distribution of the 
Sgr galaxy is dominated by a w i de peak at \Fe/H] c:^ 
-0.4 (see iMoiiaco et all l2005bl : iBonifacio et all 120061 : 
iBellazzini et al.ll2006al and references therein), extend- 
ing from [Fe/H] 1.0 to [Fe/H] ~ 0.0. Stars more 

metal poor than [Fe/H] ^ —1.2 are quite rare and, 
presumably, they would pass their core-Helium-burning 
phase as RR Lyrae or Blue Horizontal Branch stars 
([Monaco et al.l l2003h : hence Red Clump stars of Sgr 
should have [Fe/H] > —1.0, much more metal-rich than 
any M54 star. 

As a first step for the selection of two samples represen- 
tative of the Sgr,N and M54 populations, we introduce 
the photometric classification defined by the three curves 
overplotted on the CMD of Fig. [3| Proceeding from left 
(blue) to right (red) we assigned the fiag cmd in the fol- 
lowing way: cmd = 1 to the stars enclosed between the 
first and second curve, as likely M54 members, cmd = 2 
to the stars enclosed between the second and third curve, 
as likely Sgr,N members, and cmd = to stars lying to 
the blue of the first curve or to the red of the third curve. 
A check a-posteriori has shown that all the stars having 
cmd = have radial velocities incompatible with being 
members of Sgr,N or M54, thus supporting the validity 
of our photometric classification. 

As the coordinate of the centers of M54 and Sgr,N we 
adopt an = 28 3.763750° and Sq = -30.478333° from 
Novo la fc Gebh ardt (2006), and we convert to cartesian 
coordinates X,Y (in arcmin) projecting the equatorial 
coordinates of each star (a, S) on the plane of sky as in 
Ivan de Yen et all ([2006f ) 

X = -(lQ800/7r)cos{S)sin{a - ao) 

Y = (10800/7r)[sm((5)cos((5o) — cos{6)sin(6Q)cos{a — ao)] 

with X increasing toward the West and Y increasing to- 
ward the No rth. Adopting the di stance to Sgr and M54 
measured bv [Monaco et al.l ([2004( 1. 1 arcmin corresponds 
to 7.65 pc. 

3.2. Comparison between independent measures 

As there are several stars in common between the ob- 
servational sets taken with different instruments and or 
set-ups, we have the opportunity to check the consis- 
tency and the accuracy of our Vr and [Fe/H] estimates. 
In Fig. [31 we show the comparison between Vr estimates 
from the various sources: V^ are the velocities obtained 

with FLAMES/fibers, V^ are from DEIMOS/holes, and 
are from DEIMOS/slits. It can be readily appre- 
ciated that the consistency among the different sets of 
measures is excellent (i.e., AVr ~ 0.0 km/s). The typical 
accuracy, as measured from the r.m.s. of the Vr differ- 
ences is < ±2.0 km/s, that is satisfying for the present 
purpose. The actual uncertainty on the single measure 
should be a factor -\/2 smaller than the r.m.s. of AVr 
(i.e. ~ ±1.4 km/s), as the latter includes the uncertain- 
ties of both the estimates, added in quadrature. Of the 
414 velocity differences plotted in Fig. [H just ~ 10 are 
significantly larger than ±3.0 km/s: some of these may 
be associated with binary system obs erved at different 
orbital phases f see [Monaco et al]|2007[ ). 
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Fig. 3. — CMD (upper panel) and distribution in the sky (lower panel) of the 1152 stars constituting our final sample. DilTerent 
colors/symbols are associated to measures obtained with dilTerent instruments and/or set-ups. The curves in the CMD show the photometric 
selection box adopted for a first discrimination between candidate M54 or Sgr,N members. 



Fig.[5]shows the comparison between [Fe/H] estimates. 
In this case, while the consistency between measures 
from FLAMES and DEIMOS/holes spectra is very good, 
[Fe/H] from DEIMOS/slits spectra are ~ 0.08 dex larger 
- in average - than those from the other two sources. We 
were unable to find the cause of this mismatch and to 
preserve the homogeneity of the final merged dataset we 
corrected all the [Fe/H] values from DEIMOS/slits spec- 
tra by adding —0.08 to all of them. The overall accuracy 
is more than satisfying, with a typical r.m.s of 0.1 dex, 
and also in this case the outliers are rare. 

Finally Fig. [5] shows the comparison between our fi- 
nal Vr and [ F e/H] e stimates and those obtained by 
[Monaco et all (|2005bD from high-resolution FLAMES- 
UVES@VLT spectroscopy, for the seven stars in common 
between these datasets. The agreement in the radial ve- 
locity is good: if the only outlier (at AVr ^ —10 km/s) 
is excluded the r.m.s. of the difference is just 2.3 km/s 



and the average difference is zero^. Therefore the relia- 
bility and the accuracy of our radial velocity scale is fully 
confirmed. The comparison of metallicities is more dif- 
ficult to interpret and would have benefited by a larger 
number of stars in common between the two samples. 
For four of the seven stars, ranging from [Fe/H] ~ —1.5 
to [Fe/H] ~ —0.2 the agreement is more than accept- 
able, given the associated uncertainties. The other three 
stars show a considerable A[Fe/H], but one of them 
was reported t o have a tentative [Fe/H] estimate from 
[Monaco et al.l (|2005bD because its spectrum was affected 
by molecular bands making the analysis less reliable. In 
any case, none of the observed differences are so large as 
to cause us to erroneously classify a Sgr,N star as a M54 

^ Also in this case we sugges t multiplicity as the most likely 
origin of the outlier (see, again. iMonaco et al.l [20071) . Note that 
binaries are not expected to significantly affect the observed 
velocit y dispersio ns in systems like those c onsidered here (see 
Hargr eaves et anil99& : .Olszewski et al.lll996l . for detailed discus- 
sion). 
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Fig. 4. — Comparison between independent Vr es- 
timates obtained from spectra acquired with different 
telescopes/instruments/set-ups. Each panel reports the dif- 
ferences between the Vr estimates as a function of V magnitude. 
The ±3.0 km/s range around zero is enclosed by the two dotted 
lines. Upper panel: DEIMOS/holes vs. FLAMES/fibers (160 
stars); middle panel: DEIMOS/Slits vs. FLAMES/fibers (80 
stars); lower panel: DEIMOS/holes vs. DEIMOS/Slits. The 
standard deviation of each set of AVr, computed after recursive 
clipping of the very few 2.5-0- outliers, is reported in each panel. 



Fig. 6. — Comparison between the [Fe/H] and Vr es timates 
obtained in this work (t.w.) and in IMonaco et al.l Il2005bl . M05), 
for the seven stars in common. The encircled dot correspond to a 
star with molecular bands in its UVES spectrum, whose metallicity 
is classified as tentative by Monaco ot al. (2005b). In the lower 
panel, the standard deviation value reported in parentheses refers 
to the sub-sample from which the outlier at AVr — 10 km/s has 
been removed. The dotted lines encloses the range 0.0 ±0.2 dex in 
the upper panel, and 0.0 ± 3.0 km/s in the lower panel. 
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Fig. 5. — Comparison between independent [Fe/H] es- 
timates obtained from spectra acquired with different 
telescopes/instruments/set-ups. Each panel reports the differences 
between the [Fe/H] estimates as a function of V magnitude. The 
±0.2 dex range around zero is enclosed by the two dotted lines. 
Upper panel: DEIMOS/holes vs. FLAMES/fibers; middle panel: 
DEIMOS/Slits vs. FLAMES/fibers; lower panel: DEIMOS/holes 
vs. DEIMOS/Slits. The standard deviation of each set of 
A[Fe/H] is reported in each panel. 



member, or vice-versa, according to the selection crite- 
ria that are adopted below. Taking also into account 
the high degree of self-consistency among different sets 
of measures shown in Fig. [51 we conclude that also our 
metallicity scale is sufficiently reliable and accurate for 
the purposes of the present analysis. 
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Fig. 7. — Radial velocity of program stars as a function of dis- 
tance from the center of Sgr,N/M54 {upper panel), and radial ve- 
locity distribution {lower panel). The long dashed lines mark the 
range we adopted for the first selection of candidate Sgr/M54 mem- 
bers. The dotted lines enclose the (global) ±3(t range from the 
mean of the selected samples of candidates (continuous line). 

3.3. Radial velocity distribution 

The distribution of the radial velocity of all the ob- 
served stars as a function of their (projected) distance 
(r = VA^ + y^) from the center of the system is shown 
in Fig. [71 Two very different populations can be iden- 
tified in this plot: a broadly-distributed cloud of Galac- 
tic field stars showing a large dispersion around ~ 
km/s, and an abundant low-dispersion population with 
a systemic ~ 141 km/s, typical of M54 and Sgr,N. A s 
a first broad selection, and following llbata et al.l (|1997f) . 
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Fig. 8. — Zoomed view of the observed radial velocity distri- 
bution (histogram with noise error bars) around the peak of the 
Galactic field population. The thick curve is the distribution of all 
Galactic field stars with color and magnitude similar to program 
stars predicted by the R03 Galactic model in the considered direc- 
tion, shifted by —8 km/s to match the position of the observed peak 
and rescaled to minimize the value with respect to the observa- 
tions for Vr- < 100.0 km/s. The thin curve shows the distribution 
of (only) giant stars from the model, with the same normaliza- 
tion. The vertical dotted lines enclose the 100 km/s < < 180 
km/s interval. The expected number of Galactic stars falling in 
this velocity range Ncont is reported. 

we retained as possible M54/Sgr,N members all the (832) 
stars having 100 km/s < Vr < 180 km/s (dashed lines 
in Fig. [7]). We then computed the mean and the dis- 
persion of this sample: the ±3(t range around the global 
mean is enclosed within the two dotted lines. We do not 
reject the stars outside the 3ct range at this stage, since 
the velocity distribution is different at different radii (see 
below), and so a more refined choice is to exclude > 3a 
outliers in any considered radial bin. However, in Sect. 4 
we will see that the adopted bin-by-bin 3ct rejection cri- 
terium will exclude all of these stars. 

Fig. [5] shows that we expect some (limited) contamina- 
tion from Milky Way stars even in our Vr selected sample. 
The thick continuous line superposed on the observed 
histogram of the radial velocities zoomed on the Galac- 
tic component, shows the excellent agreem ent between 
the ob servations and the predictions of the I Robin et al.l 
(|2003D Galactic model, once it has been shifted by —8.0 
km/s to match the observed peak and rescaled to mini- 
mize for Vr < 100.0 km/s. The model suggests that 
up to ~ 30 Milky Way stars may be present even in the 
relatively narrow Vr window we have adopted to select 
M54/Sgr,N stars. It is interesting to note that among the 
stars having 100 km/s < Vr < 180 km/s but lying outside 
of the {Vr)±3a range, seven have Vr < {Vr) — 3a and four 
have Vr > (Vr) + 3(T, in agreement with the trend with 
Vr predicted by the model. Fig. [8] shows also that most 
of the Galactic interlopers are exp ected to be giants, ac- 
cording the the lRobin et al.l ()2003f ) Galactic model. How- 
ever, it is likely that such stars should be - in general - 
much closer to us than Sgr: the model finds that the con- 
taminating giants have a mean distance of 8.7 kpc (corre- 
sponding to A{M—m) = +2A mag, with respect to Sgr), 
90% of them have D < 11 kpc and 99% of them have 



D < lA kpc, i.e. much lower than Dsgr — 26.3 kpc. For 
this reason the difference between their V magnitude and 
Vhb of Sgr must be a bad proxy for their gravity. The 
typical spectroscopic metallicity error incurred by over- 
estimating the distance modulus of these Galactic stars 
by 2.4 mag is A[Fe/H] = —0.58, so the stars will appear 
to be less metal-rich than they are in reality while their 
color and magnitude mimic those of genuine Sgr/M54 
stars. As a consequence, for a Galactic star the metallic- 
ity obtained from the Calcium Triplet (CaT) technique 
would normally be much different from what would be 
expected from the magnitude and color measurements. 
Therefore, there is some hope to exclude these Galactic 
interlopers from our final samples, identifying them by 
their "odd" color-metallicity combination, as we will do 
in our final selection (Sect. 3.5, below). As the mean 
metallicity of contaminating giants is {[Fe/H]) — —0.75 
and 90% of them have [Fe/H] < —0.33 it is quite likely 
that many of them would have a measured metallicity 
lower than the lower threshold we will assume for Sgr,N 
members ([Fe/H]> —0.8). On the other hand those spu- 
riously showing metallicities compatible with M54 will 
probably be too red to be selected as possible cluster 
members (Sect. 3.5 and Fig. [3]). 

Finally it is worth noting that the spectra of all the 843 
stars that passed the selection in Vr have Signal to Noise 
ratio (per pixel) S/N > 12; 300 of them have S/N > 50 
and 63 of them have S/N > 100 (see Fig. Oi). 

3.4. Metallicity distribution 

The metallicity distribution of the velocity selected 
stars is shown if Fig. [5^. The presence of two popula- 
tions is very clear: there is a metal poor peak around 
\Fe/H] :^ -1.5 that must be domina ted by M54 stars 
(jBrown. Wallerstein fc Gonzaled Il999l ) and a broader 
distribution extending from [Fe/H] > —1.0 to super- 
solar metal hcity corresponding t o the main population of 
Sgr.N (seelMonaco et all l2005bl: iBellazzini etani2006at 
iBonifacio et al.ll2006l and references therein). As men- 
tioned earlier, in the present context we will use the 
metallicity just as a further means to select samples of 
M54 and Sgr,N members that are as clean as possible, so 
as to allow further analysis of their kinematic properties. 
Therefore, we limit the discussion of the metallicity to 
the brief considerations listed below. 

1. We confirm that the bulk of Sgr stars have 
[Fe/H] > —1.0, belonging to a broad distribu- 
tion peaking around [Fe/H ] = —OA an d reach- 
ing super-solar metalliciti es (iMonaco et a l. 2005^; 
IBonifacio et~all [2006( ) . iMonaco et all (|2603i ) esti- 
mated that in the main body of the Sgr galaxy old 
metal poor stars {[Fe/H] < —1.2) should provide 
less than 12% of the whole stellar content. Sev- 
eral arguments suggest that this fraction should 
be even lower in the very central region consid- 
ered here (see lAbIdll200ll : IBellazzini et aTll2006bl : 
ISiegel et al.l l2007l and references therein). How- 
ever, we note that there are four stars having 
[Fe/H] < —2.0 in our sample (see Fig. ^ that are 
probably too metal poor to be members of M54 
and might be part of the metal-poor population of 
Sgr. 
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Fig. 9. — Selection of clean samples of M54 and Sgr,N members. Panel a: metallicity distribution of the 843 stars that survived 
the selection in radial velocity. The superposed curve is the sum of two Gaussian distributions having mean and standard deviation 
{{[Fe/ H]) , a) = (—0.45,0.28) and (—1.45,0.17), representing the main populations of Sgr,N and M54, respectively. The dotted histogram 
is the observed distribution for stars having R < 3.0'. Panel b: metallicity vs. I magnitude; open triangles are stars having cmd = 2, 
i.e. photometrically selected as likely Sgr members (see Fig. [Sj; open pentagons are stars with cmd = 1, i.e. photometrically selected as 
likely M54 members. The horizontal lines enclose the Sgr RC magnitude range (16.6 < I < 17.3). The vertical lines mark the adopted 
thresholds in metallicity: —1.8 < [Fe/H] < —1.1, for the M54 sample; —0.8 < Fe/H] < +0.2 for the Sgr,N sample. Panel c: metallicity 
vs. distance from the center of M54. The symbols are the same as in the above panel, except for asterisks superposed on open pentagons: 
these are stars that lie on the RGB of M54 in the CMD but in fact belong to the superposed RC of Sgr, i.e. they have 16.6 < I < 17.3 
and a metallicity typical of Sgr members. Panel d: Signal to Noise (per pixel) of the spectra as a function of I magnitude. All the spectra 
have S/N > 12. 



The peak of the metalhcity distribution corre- 
sponding to M54 is best-fitted by a Gaussian 
curve having mean {[Fe/H]) — —1.45. This is 
sh ghtly higher -1-0.1) than w h at wa s found 
by Brown. Wa llerstein fc Gonz ale3 ()1999D and by 



iDa Costa fc A rmandroff (199 5), w hile it is in good 
agreement with lArmandrofl (119891 ) . Note that any 
possible small shift in the zero point of our metallic- 
ity scale does not affect the accuracy of the metal- 
licity ranking, that is ~ ±0.1 and it is the relevant 
figure in the present context. 

3. Based on the analysis of th e Color - Magnitude 
distribution of RGB stars, ISaraiedini &: LavdenI 



(|1995D concluded that there is an intrinsic metal- 
licity spread among M54 stars of aint{[Fe/ H]) = 
0.16 dex. Using spectroscopic metallicities ob- 
tained wi th the CaT technique f or fiv e M54 
members, iDa Costa fc Armandrofj ()1995[ ) found 
further support for this hypothes is (but see 
iBrown. Wallerstein fc Gonzalej [19991 ). Our much 
larger sample of CaT metallicities provides new 
supporting evidence for this possibility. The sam- 
ple with -1.0 < [Fe/H] < -2.0 has an ob- 
served standard deviation of 0.17 dex, larger than 
what is expected from the statistical scatter (see 
Fig. [51 above). Deconvolving the (internal) r.m.s. 
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Fig. 10. — Upper left panel: CMD of stars selected to belong to t he final Sgr,N (open triangles) and M54 (open pentagons) samples, 
superposed on the CMD of all stars having R < 9' from lMonaco et al.l (2002, small dots) . The thin horizontal lines enclose the Red Clump 
of Sgr (16.6 < I < 17.3). The position in the sky of the selected stars relative to the center of Sgr,N/M54 is plotted in the lower left panel, 
the symbols are the same as above. The radial distributions of stars selected on the Sgr,N sequences (continuous histogram) and on the 
M54 sequences (dotted line) are plotted in the lower right panel The upper right panel shows the uncertainties of the Vr estimates as a 
function of I magnitudes for the selected stars. 



scatter (0.11 dex, see Fig. [SJ from the observed 
width of the distribution we obtain an estimate of 
the intrinsic spread of aint{[Fe/H]) = 0.14 dex, 
in good agreement with previous estimates. It 
seems very unlikely that the width of the M54 
peak is significantly contaminated by Sgr stars, as 
the width remains unchanged if we limit the sam- 
ple to the innermost 3' where M54 should dom- 
inate the population mix. If finally confirmed 
by high resolution spectroscopy of a large sample 
of stars, the purported metallicity spread would 
constitute another element of similarity between 
M54 and other large-size very bright clusters like 
uj Cen, Gl and the like, that are also suspected 
of being remnan t nuclei of disrupted galaxies 



(see Mackev & van den Bcrgh 2005; Fcdcrici et al 
12007, : ■Catelan..2007. : Maschenko fc Sills 2005, : iBoked 



|2007| . and references therein). 

We think that the presently available metallicity esti- 
mates provide a metallicity ranking that is more than 
sufficient for the purposes of the present paper, i.e. (a) 
to discriminate M54 stars from Sgr,N stars in the mag- 
nitude range where the two population overlaps in the 
CMD (16.8 < I < 17.3), and (b) to exclude possible in- 
terlopers as stars having the "wrong" metallicity for their 
color. The details of the adopted selections are described 
in the following sub-section. 

3.5. The final selections 

The scheme of our finally adopted selections is illus- 
trated in Fig. inbjC. Our aim is to have a sample of M54 
stars and a sample of Sgr,N stars as clean as possible 
from any kind of interlopers. The adopted criteria are 
the following: 
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• We accept as members of the M54 sample the stars 
having cmd = 1 and -1.8 < [Fe/H] < -1.1. 

• We accept as members of the Sgr,N sample the 
stars having cmd = 2 and -0.8 < [Fe/H] < +0.2, 
or having cmd = 1 and -0.8 < [Fe/H] < +0.2 if 
16.6 < I < 17.3 (Sgr RC stars superposed on the 
RGB of M54). 

The metallicity separation avoids mixing between sam- 
ples. The requirement that the cmd flag and the metal- 
licity broadly agree in assigning the membership is likely 
to exclude most of the Galactic interlopers. In the region 
in which the RC of Sgr,N and the RGB of M54 overlap, 
the cmd flag does not provide any discrimination and 
we rely just on the metallicity. We have adopted quite 
conservative criteria to obtain very reliable samples of 
genuine Sgr,N and M54 stars. As a further check, all the 
results presented in the following about the kinematics of 
Sgr,N have been verified to hold also when subsamples in 
which stars in the RC region were excluded are adopted. 
The main properties of the selected samples are shown 
in Fig. ]M 

With the selection described above we select a sample 
of 425 very likely M54 members and of 321 very likely 
Sgr,N members. In the process of clipping 3cr velocity 
outliers in any considered radial bin (which will be per- 
formed in Sect. 4, below), we will further exclude eight 
stars from the M54 sample and three stars from the Sgr,N 
sample. Hence the actual detailed analysis of the kine- 
matics of the two systems will be performed on 417 stars 
for M54 and on 318 stars for Sgr,N. 

4. THE KINEMATICS OF M54 AND SGR.N: 
OBSERVATIONAL FACTS 

In Tab. |4] we list the estimates of some relevant physi- 
cal parameters of M54 that are available in the literature. 
For the structural parameter s (r,., rt, etc.) we will pref- 
erenti ally adopt those from iTrager. King fc Diorgovskil 
since the best-fit model proposed by these 
authors is in good agreement with the profile obtained 
from star counts in the outer regions of the cluster 
by M05a. However in some cases we will consider the 
effects of the adoption of different sets of parameters 
on our results. It is reassuring that most authors 
derived quite similar parameters from very different 
datasets. All the available mass estimates for M54, 
also obtained with different methods and from different 
datasets, ran g e betw ee n 1.0 x lO^Mp and 2.0 x lO^Mfr, 
JlUingworth' '1976'; 'Mandushev, Staneva & Spasova' 
1P91; Prvor & Mevlan 1993; Mackcy & Gilmorc" 2003; 
McLaughl in fc van der Marelll2005 ). 

In Tab. [S] we report the average kinematical properties 
of the M54 and Sgr,N samples as a whole. The errors 
on all the reported quantities have been obtaine d by a 
jackknife bootstrapping technique ()Luptonlll993f ). The 
plain average Vr of the two samples differs by 1.5 ± 0.7 
km/s, that is, just a ~ 2cr difference. However an aver- 
age with iterative 2-a clipping seems a much more ap- 
propriate estimator of the systemic velocity since it is 
very robust to outliers and is more effective in estimat- 
ing the mode of the samples. The sigma clipping process 
converges very fast in both cases, leaving very rich fi- 
nal samples {N2a = 246 and 348, for Sgr,N and M54, 
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Fig. 11. — Velocity dispersion profile of M54 stars. The upper 
panel shows the Vr distribution as a function of distance from the 
cluster center for individual stars of the M54 sample. Only stars 
plotted as dots encircled by open pentagons are retained for the 
computation of a in the various radial bins: small dots alone are 
stars rejected only because they are "local" 3cr outliers of the bins, 
crosses are stars that would have been rejected also as Scr outliers 
of the whole Sgr,N + M54 sample (the global ±3(t range is enclosed 
by the long-dashed lines). The vertical lines display the adopted 
independent bins, of variable size. The global mean is marked by 
the continuous horizontal thick line. The lower panel displays the 
actual Velocity dispersion profile. The large filled pentagons are 
the dispersions estimated in the corresponding bins displayed in the 
upper panel, with their bootstrapped errors. The number of stars 
per bin is also reported below the points. The small filled pentagons 
are the estimates in the additional, partially overlapping, bins. The 
op en pentagon is the estimate of a at the center of M54 obtained 
by llUingworthl 1)19761 ') from integrated spectroscopy. 

respectively). The difference in the clipped mean veloc- 
ities is a mere 0.8 km/s. We conclude that the Sgr,N 
and M54 do coincide also in line of sight (l.o.s.) velocity, 
a further proof in support of their spatial coincidence. 
Many authors have noted before that the systemic ve- 
locity of Sgr and M54 are very similar (llbata et al.lll997t 
iDa Costa fc Arm androff '1995': 'Mona co et al.ll2005b[ ): we 
think we have conclusively established here that M54 and 
the nucleus of Sgr have the same systemic l.o.s velocity, 
within ~ ±1.0 km/s. 

The velocity dispersions of the whole Sgr,N and M54 
do not seem very different, at a first glance. However, it 
should be recalled here that in the radial range covered 
by our data (r < 9'), the surface brightness of M54 and 
Sgr,N declines by more than ^ 8 mag/arcsec^. Hence, in 
the present context, what is really relevant is the com- 
parison between velocity dispersion profiles as a function 
of distance from the center of the system, which will be 
presented in the following section. 

4.1. Velocity dispersion profiles 

In the upper panel of Fig. [TT] we show the distribution 
of the M54 stars in the V,- vs. r plane. We have divided 
the sample into two sets of six and five independent ra- 
dial bins of different size, respectively, in order to keep 
the number of stars per bin as high as possible while 
maintaining the highest degree of spatial resolution (the 
primary set of independent bins corresponds to the odd 
rows of Tab.[6l the secondary set corresponds to the even 
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TABLE 4 
Literature data for M54. 











Vt 


\Fc /HI 

[rC/ JTLJ 


method 


Ref. 


arcmin 


arcmin 




mag/arcsec^ 


mag 


dex 






0.11 


7.4 


1.83 


14.90 


7.61 
7.68 




g.b. ap. phot, 
g.b. ap. phot. 


lUinffworth & IlUnsworth (1976") 
Peterson (1993) 


0.11 


7.4 


1.84 


14.75 




-1.55 
-1.55 


g.b. ap. phot. 
CaT spectr. 
HR spectr. 


Traser. Kine & Dioreovski (1995) 
Da Costa & Armandroff (1995) 
Brown, Wallerstein & Gonzalez (1999) 


<0.11 
"o.09 


9.9 


2.04 


14.35 






HST phot, 
g.b. ap. phot. 


Mackev & Gilmore (2003) 
McLaughUn & van der Marcl ( 2005) 


0.05 






14.12 






HST phot. 


Noyola & Gebhardt (2006) 



Note. — g.b. ap. phot. = ground based aperture photometry; HST phot. = integrated photometry and/or star counts from Hubble 
Space Telescope data; CaT spectr. = medium resolution spectroscopy using the infrared Calciu m Triplet as metallicity indicator; HR 
spectr . = High Resolution spectroscopy and elemental abundance analysis. See also iWebbinh l ll985ii ). rHarri£l Ill996i) and , Mackev fc Gilmor^ 
1120031) for collections of literature data on M54. 



TABLE 5 

Average kinematic properties of the M54 and Sgr,N 

SAMPLES. 







{Vr) 


a 


a* 


N2a 








km/s 


km/s 


km/s 




km/s 


Sgr,N 


321 


139.4 ±0.6 


10.0 ±0.5 


9.6 ±0.4 


246 


139.9 


M54 


425 


140.9 ±0.4 


9.3 ±0.5 


8.3 ±0.3 


348 


140.7 



Note. — N2o- = number of stars in the sample once the iterative 
2(T clipping algorithm to compute {Vr)2<7 has converged, i.e. there 
is no more 2cr outlier. Note that at each iteration the algorithm 
adopts the current value of 2cr that is, in general, significantly 
lower than the global value reported in the table, a* are the global 
values of the dispersion after the exclusion of the local > 3(t outliers 
performed in Sect. 4.1. They are computed from the final cleaned 
samples of 318 (Sgr,N) and 417 (M54) stars, as already reported 
in Sect. 3.5, above. 

rows of the table) . The primary bins are enclosed by the 
vertical lines. In each bin we computed the average Vr 
and the velocity dispersion a, with their uncertainties 
(derived with the "jackknife" method, as above). An it- 
erative 3ct clipping algorithm was applied bin by bin. As 
we proceeded from the innermost bin to the outer ones, 
any star rejected in a given bin by the clipping algorithm 
was excluded from the sample. The eight rejected stars 
are clearly indicated in the plot. 

The derived velocity dispersion profile is reported in 
the lower panel of Fig. [TT] (large filled pentagons, primary 
set, small filled pentagons, secondary set) and in Tab-IH 
The table reports also an alternative estimate of a ob- 
tained with the Gaussian maximum-likelihood method 
described bv IWalker et al.l 12006) . essen tially equivalent 
to that adopted bv iMartin et al.l (|2007D . It is remark- 
able that in all cases the estimates obtained with the two 
methods differ by much less than the reported uncertain- 
ties (in all cases by < 0.3 km/s). The profile is comple- 
mented with the central estimate obtained by lUingworthl 
(|1976l ) from integrated spectroscopy of the cluster core''. 
The profile shows a steep decrease from a = 14.2 km/s 
at the center to a ~ 5.0 km/s at r = 3.5'. Then it begins 
to grow gently up to cr ~ 9 km/s in the last bin (r ~ 7'). 

Fig. [12] shows the velocity distribution and dispersion 
profile for the Sgr,N sample, with the same arrangement 
as in Fig. [Til above. The absence of any obvious trend 

Here we adopt the estimate of central velocity dispersion 
(TO = 14.2 km/s obtained by lUingworth by correcting the value 
he measured from integrated long-slit spectroscopy of the cluster 
center. 



TABLE 6 

Velocity dispersion profile for M54. 



ri 




rm 


a 


eo- 


0"mi 


eo- 


N 


arcmin 


arcmin 


arcmin 


km/s 


km/s 


km/s 


km/s 




0.1 


1.0 


0.68 


10.2 


0.9 


10.1 


1.0 


95 


0.5 


1.5 


0.96 


10.2 


0.7 


10.1 


1.1 


131 


1.0 


2.0 


1.50 


8.6 


0.5 


8.5 


1.2 


107 


1.5 


2.5 


1.98 


7.4 


0.4 


7.2 


1.4 


101 


2.0 


3.0 


2.47 


7.4 


0.5 


7.3 


1.3 


86 


2.5 


3.5 


2.95 


6.2 


0.6 


6.1 


1.3 


66 


3.0 


4.0 


3.38 


5.3 


0.7 


5.1 


1.2 


42 


3.5 


5.0 


4.29 


6.3 


0.5 


6.2 


1.1 


49 


4.0 


6.0 


4.83 


6.7 


0.5 


6.6 


1.1 


53 


5.0 


7.0 


5.89 


7.0 


0.9 


6.9 


0.9 


36 


6.0 


9.0 


6.87 


8.8 


1.5 


8.6 


0.7 


34 



Note. — ri and rf are the limits of the bins, rm is the radius 
of the middle of the bins. N is the number of stars in the bin, 
after the rejection of 3(7 outliers, a^i is the velocity dispersion 
estimated with a G aussian maximum likelihood method as done in 
IWalker et all 1 1200^) . 



TABLE 7 

Velocity dispersion profile for Sgr,N. 



ri 




rm 


(7 


eo- 






N 


arcmin 


arcmin 


arcmin 


km/s 


km/s 


km/s 


km/s 




0.8 


2.0 


1.41 


9.1 


0.9 


8.9 


0.7 


41 


2.0 


3.0 


2.55 


8.8 


1.2 


8.6 


0.7 


38 


3.0 


4.0 


3.54 


10.2 


0.9 


10.0 


0.7 


56 


4.0 


5.0 


4.47 


8.8 


1.0 


8.5 


0.9 


43 


5.0 


6.0 


5.46 


9.6 


1.1 


9.3 


0.7 


40 


6.0 


7.0 


6.52 


8.7 


0.9 


8.5 


0.7 


43 


7.0 


9.0 


7.69 


10.3 


1.1 


10.1 


0.7 


56 



Note. — and rf are the limits of the bins, rm is the radius 
of the middle of the bins. N is the number of stars in the bin, 
after the rejection of 3cr outliers. (t„; is the velocity dispersion 
estimated with a G aussian maximum likelihood method as done in 
IWalker et all 1(2003 ). 



of velocity dispersion with radius allowed us to adopt a 
single set of bins of nearly uniform size. The dispersion 
in each bin has been estimated exactly as in the case of 
M54, described above. The two innermost stars of the 
sample (crosses) have been excluded from the analysis as 
they lie in a too poorly sampled radial region (r < 0.8') 
and may affect the computation of the mean radius of 
the first bin. However, the inclusion of these stars in the 
innermost bin changes the velocity dispersion estimate 
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Fig. 12. — The same as Fig. Illl for the Sgr,N sample (triangles 
replaces pentagons here). The two innermost point are rejected for 
reasons of uniformity of the sample. In the lower panel we also 
report the velocity dispersion of the whole sample as a continuous 
line, ± the associated error (dotted lines). As u does not appear 
to change with radius we use only the reported seven independent 
bins. 

by ~ 0.2 km/s. 

It is quite obvious from the inspection of Fig. [T^that 
the velocity dispersion profile of Sgr,N is completely flat 
over the whole radial range explored by our data. The 
lower panel shows that the velocity dispersion of each 
radial bin is in good agreement with the velocity disper- 
sion of the whole sample after the clipping of 3fT outliers, 
cr = 9.6 ± 0.4 km/s. The maximum observed difference 
among bins is 1.6 ± 1.4 km/s; the maximum difference 
between the dispersion in a given bin and the dispersion 
of the whole sample is —0.9 ± 1.0. Possibly, there may 
even be a weak tendency toward lower velocity disper- 
sions in the inner regions: in the range 0.0' < r < 5.0' 
wc find cr = 9.2 ±0.5 km/s (180 stars), while in the range 
5.0' < r < 9.0' we find a = 10.1 ± 0.6 km/s (138 stars); 
however th is difference is clearly not statistically signif- 
icant (see iGeha. Guhathakurta fc van der Marell [200l : 
iValluri et al.l l2005l for some example of the varietv of 
velocity dispersion profiles observed in dE,N). 

It is rather clear that the nucleus of Sgr and the M54 
cluster have different velocity dispersion profiles. M54 
displays the typical behavior of a globular cluster, i.e. the 
system becomes increasingly kinematically hot toward its 
central regions while Sgr,N has a uniformly flat profile, 
though the surface brightness profiles of the two systems 
are broadly similar. Fig. [T3] shows the comparison be- 
tween the profiles in deeper detail. The inner growth of 
cr of M54 stars makes the two profiles overlap at r ~ 1.5'; 
the innermost two points of the M54 profile are hotter 
than any Sgr,N point. The outermost point of the cluster 
profile does match the dispersion of Sgr,N: the velocity 
dispersion estimate here comes from the least populated 
bin as well as the most (possibly) prone to contamination 
by old and metal-poor stars belonging to Sgr (see below 
for a more detailed discussion). These are the reasons 
why the comparison between the two samples taken as a 
whole does not reveal any striking difference. However, 
in the intermediate radial range the M54 population is 



Fig. 13. — Upper panel: Velocity distribution of the M54 and 
Sgr,N samples in the radial range in which the difference in their 
(7 reaches the maximum statistical significance, according to the 
F test (see Tab. |9] Lower panel: Comparison of the velocity dis- 
persion profiles of M54 (continuous line and pentagons) and Sgr,N 
(dotted line and triangles). 



significantly colder than that of Sgr,N and the statistical 
significance of this result is very strong. 

Tab. H] reports the results of non-parametric 
Kolmogorov-Smirnov tests performed on different radial 
ranges. In all the considered ranges the probability that 
the two samples are drawn from the same parent popu- 
lation is smaller than 2% and it is smaller than 0.2% in 
the 1.5' < r < 6.5' range. In the present case, however, 
the F statistic is the most appropriate mean to com- 
pare statistically the two samples {F = o-g^^/a^j^^, see 
lBrandtl[l970l . and references therein), as the correspond- 
ing F test evaluates the probability that two samples of 
given F value are extracted from Gaussian distributions 
having the same a. Tab. [9] shows that in our case this 
probability is lower than 0.3% in all the considered radial 
range and it is as small as < 0.01% in several intermedi- 
ate ranges, including 1.5' < r < 6.5'. While such a large 
quantitative difference may not be so obviously apparent 
from the dispersion profile, it is clear at a first glance 
when looking at the comparisons in the X - and in 
the K - Vr planes shown in Fig. [Ml 

Hence, the direct comparison of the velocity distribu- 
tions of the two samples clearly establishes that the two 
populations have very different kinematic properties, in- 
dicating that the motion of M54 and Sgr, N stars is driven 
by different gravitational equilibria. 

Moreover, iGilmore et al.. (2007D clearly states that a 
velocity dispersion profile monotonically declining from 
the center outward "...is an unavoidable requirement for 
any mass - follows - light system. . . ". In the present 
case, this means that while M54 behaves like an ordinary 
mass - follows - light "purely baryonic" self-gravitating 
star cluster (at least in its innermost region, containing 
most of its light/mass; see also Fig. 14 and Sect. 4.2, be- 
low), the flat velocity dispersion profile of Sgr,N implies 
that the mass distribution driving the kinematics of the 
nucleus is significantly different from the distribution of 
the stars. In principle, a radial gradient in the velocity 
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Fig. 14. — Direct comparison of the velocity and spatial distribu- 
tions of M54 and Sgr,N. Upper panels: Vr as a function of X for the 
M54 (left panel) and for the Sgr,N (right panel) samples. Lower 
panels: Vr as a function of R for the same samples. The curves in 
the lower left panel are the ili3fT profiles t he best-fit King model 
for M54 IjTrager. King &: Diorgovskilll995i ). approximately repre- 
senting the envelope of the allowed velocities for bound members 
of M54. 



TABLE 8 

KS TEST. Probability that the Sgr.N and the M54 samples 

ARE DRAWN FROM THE SAME PARENT POPULATION. 







Nm54 




Dks 


P% 


0.0 


9.0 


417 


318 


1.649 


0.87 


0.8 


8.0 


353 


300 


1.531 


1.84 


2.0 


7.0 


199 


219 


1.624 


1.02 


2.0 


6.0 


181 


177 


1.698 


0.63 


1.5 


6.5 


247 


213 


1.817 


0.13 



Note. — KS = Kolmogorov-Smirnov. and rf are the limits 
of the considered radial ranges, N^ampic is the number of stars of 
the sample in the considered radial range. Scr outliers have been 
excluded from the samples. 



TABLE 9 

F TEST. Probability that the Sgr,N and the M54 samples 

ARE DRAWN FROM GAUSSIAN DISTRIBUTIONS HAVING THE SAME cr. 





rf 


Nm54 




F 


P% 


0.0 


9.0 


417 


318 


1.345 


0.23 


0.8 


8.0 


353 


300 


1.448 


0.04 


2.0 


7.0 


199 


219 


2.010 


<0.01 


2.0 


6.0 


181 


177 


2.019 


<0.01 


1.5 


6.5 


247 


213 


1.978 


<0.01 


2.5 


3.5 


65 


51 


2.546 


0.02 


1.5 


3.5 


166 


80 


1.862 


0.04 



Note. — and r-f are the limits of the considered radial ranges, 
^sample number of stars of the sample in the considered 

radial range. Scr outliers have been excluded from the samples. F 
is computed as F = crl^^ j^/alj^^. 




Sgr galaxy, that is found to be flat out to large radii 



(llbata et al.1 1997D . as typical of dwarf spheroidals (see 



Fig. 15. — Velocity dispersion profile of the innermost ~ 100' of 
the Sgr galaxy. Filled triangles are the r < 10' estimates obtained 
here, the skeletal symbol is from the M05a stars having [Fe/H] < 
— 1.0 (to avoid possible M54 members), open triangles are from 
Table 2b of Ibata ct aO 1(19971 ) (CTIO-ARGUS radial velocities). 
The horizontal linos are the same as in Fig. 1121 Note that this 
composite profile covers less than half of the core radius of the Sgr 
galaxy (re ~ 220' iMaiewski et al.ll2003D . 

anisotropy affecting only one of the two systems may be 
invoked to explane the observed difference in the veloc- 
ity dispersion profiles. As this would imply a correlation 
of the anisotropy variation with metallicity, we regard 
this hypotesis as rather "ad-hoc" and we don't discuss it 
anymore, in the following. We will re-consider the case 
in more detail in a future contribution (Ibata et al., in 
preparation) . 

Finally, the velocity dispersion profile of Sgr,N does 
not differ significantly from the overall profile of the 



IWalker et al.ll2006t iMufioz. Maiewski fc JohnstonI [20071 
and references therein). In Fig. [T3] we have combined 
the inner 0.8' < r < 9' profile obtained here, with esti- 
mates of the dispersion at different radii out to r ~ 100' 
obtained fro m publicly availa ble data from the litera- 
ture (M05a. llbata et al.lll997f) in the same way as done 
for M54 and Sgr,N, above. Note that r ~ 100' is less 
than half of the core radius of the whole Sgr galaxy 
(|Maiewski et al.l 120031 ). All the observations are consis- 
tent with a constant a in the considered range and, in 
particular, there is no apparent transition in the profile 
at the onset of the nuclear overdensity of stars (r < 10'). 
This fact strongly suggests that the inner kinematics of 
Sgr, including its nucleus, are dominated by the poten- 
tial set by a distribution of unseen (Dark) matter, not by 
M54 as in the scenario (b) described in Sect. 1., above. 

4.2. Comparison with models 

In Fig. [in] the observed velocity dispersion profile of 
M54 is compa red wi t h the theoretical profiles of single- 
mass isotropic iKingl (|1966f ) models that have been pro- 
posed by various authors as best-fits to the surface 
brightness profil e of the cluster (see Tab. |4| . The mod- 
els proposed by 'Trager. Kin g fc Diorgovskil (Il995|, con- 
tinuous line) and bv [McLaughlin fc van der Marell ll2005l . 
dashed line) provide a reasonable description of the data 
out to r < 4', while there is a clear tendency of the data 
to be hotter than the models for r > 4'. The outermost 
observed point, at r ~ 7' < r^, displays a velocity dis- 
persion significantly higher than what is predicted even 
by the hottest model [dofes ^ o-rm] ^-,1 — 5.8 ± 2.0 km/s. 



16 



Bellazzini et al. 




Fig. 16. — Comparison of the observed velocity dispersion pro- 
file of M54 with King models that best fit the Surface Brightness 
profile of the cluster, as proposed by various authors (see Tab.|4]l, 
and normalized to ctq = 14.2 km/s. Continuous line: Tc = 0.11' 
and C = 1.85, from Tragcr, King & Diorgovski (1995). Dotted 
line: rc = 0.053' from Novola k, Gcbhardt (,200a) and C = 2.10, 
obtained by adopting the estimate of the tidal radius from Trager 
et al., since it is not provided by Noyola &: Gebhardt. Dashed 
line: = 0.09' and C = 2.05, from [McLaughlin &: van der Marell 
ll2005h . 

There are two possible explanations for this unexpected 
rise of the velocity dispersion curve near the tidal radius 
of the cluster. 

First, while our selection criteria prevents any contam- 
ination of the Sgr,N sample by M54 stars, since cluster 
stars with [Fe/H] > —0.8 simply do not exist, it is well 
known that the Sgr galaxy d oes host a minority of stars 
with [Fe/H] < -1.0 (<12%. [Monaco et al.l[2003l ). hence 
some degree of contamination of the M54 sample by Sgr 
stars is expected (see Sect. 3.1) and the less populated 
outer bins of the M54 profile are clearly the most sensi- 
tive to the effect of the inclusion of contaminants. Rough 
estimates indicate that the contamination by metal poor 
Sgr stars should amount to less than 3% in the inner- 
most 2', less than 10% for 2' < r < 4', but it can raise 
to < 20% for 4' < r < 6' and be even larger for r > 6'. 
So, contamination by metal-poor Sgr stars seems a viable 
explanation for the raise of the dispersion profile of M54 
such that it becomes similar to that of the Sgr,N sample 
in the outermost, and presumably most contaminated, 
bin. 

Alternatively, it may be conceived that, during its spi- 
raling toward the center of Sgr,N (see Sect. 5.2, below), 
M54 suffered some tidal harassment from the host galaxy 
(Sgr). In this case, some stars in the last bins would 
have been stripped from the outskirts of the cluster and 
heat ed up to the same dispersion of the su rrounding field 
(see iMunoz. Maiewski &: Johnstonll2007l and references 
therein). In the lower left panel of Fig. [14] the ±3(7 dis- 
persion profiles of the best-fittin g King model for M54 
([Trager. King &: Diorgovskilll995D are superposed on the 
M54 sample in the R vs. Vr plane: the good match be- 
tween the shape of the bulk of the observed distribution 
and the model profile out to the tidal radius strongly 
suggests that most of the velocity outliers at large radii 



Fig. 17. — Comparison of the observed velocity dispersion profile 
of Sgr,N with theoretical models. Upper panel: comparison with 
King models of rc = 0.05' (best fit value) and various C (thin 
continuous lines labeled with their C value), normalized at the 
central velocity dispersion of M54 (co = 14.2 km/s). The C ~ 1.9 
case corresponds to the model that best-fits the observed SB profile 
and is plotted with a thicker line. The dotted curves show the effect 
on the profile of the assumption of Tc and rt values much larger 
than the best fit value (of a factor x2). Lower panel: the observed 
velocity dispersion profile is compared with the profile of the NFW 
model adopted in the N-body simulation cll5_b (see Tab. llO|l . The 
profiles at the beginning (continuous line) and at the end (dotted 
line) of the simulation are shown. 

are likely not gravitationally bound to the cluster, inde- 
pendent of their origin (unrelated metal poor Sgr stars 
or former M54 members that have been tidally stripped). 

While we consider the "contamination" hypothesis as 
more likely, the "tidal" hypothesis is a very fascinating 
possibility and it cannot be dismissed with the data we 
have presently in hand (see also Sect. 4.3, below). 

Turning back to the inner regions of the cluster, the 
main conclusion we can draw here is that the observed 
velocity dispersion profile is consistent with the expected 
kinematics of the King models that best fit the light dis- 
tribution of the cluster. Hence, M54 has the kinematics 
of an ordinary self- gravitating globular cluster, at least 
in its innermost ~ 25 pc, enclosing more than 90% of the 
cluster light/mass. 

4.2.1. The case of Sgr,N 

In Fig. [T7| the observed velocity dispersion profile of 
Sgr,N is compared with various theoretical models. In 
the upper panel we compare King's models with rc — 
0.05' and various values of the concentration parameter 
C (the best-fit model for the SB profile corresponds to 
C=1.9). All the theoretical profiles have been normalized 
at the central velocity dispersion of M54: this is a quite 
arbitrary choice, the family of profiles can be shifted up 
and down according to the preference of the reader. It 
is clear that, independent of the adopted normalization, 
none of the plotted profiles provides a satisfactory match 
to the observations. The dotted profiles, reported here 
for comparison, show the effects of factor 2x changes 
in the adopted values of rc and rt- Even if such large 
values were compatible with the observations presented 
in Sect. 2, above (and they are not), the corresponding 
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profiles remain unable to reproduce the observed disper- 
sions. 

One could be tempted to suggest that any King model 
with sufficiently large scale to fit the whole surface pro- 
file of the Sgr ga laxy (as the C=0.90, = 224' model by 
iMajewski et ali2003) would have a nearly flat dispersion 
profile on the relatively small radial range studied here 
(r < 9'), thus providing a reasonable fit to the observed 
profile. This is true, but these models would completely 
fail in reproducing the observed Surface Brightness pro- 
file in the same range by more than 10 mag/arcsec^ (see 
Fig. [21 above), since they lack any nuclear overdensity at 
their center. 

Hence, contrary to the case of M54, King models that 
fit the surface brightness profile of Sgr,N appear unable 
to reproduce its kinematics: this further supports the 
conclusion that M54 and Sgr,N are systems of very dif- 
ferent nature*. In particular, as already noted, mass - 
follows -light models are incompati ble with the observe d 
velocity dispersion profile of Sgr,N l|Gilmore et al.ll2007t ). 

In the lower panel of Fig.[T7]the observed profile is com- 
pared to those obtained f r om th e N-body realization of a 
ikavarro. Frenk fc Wh"it3 ()1996l . hereafter NFW) model 
of the suite that is described in Sect. 5, below (model 
NFW3, Tab. fTO]). The continuous line is the inner veloc- 
ity dispersion profile of the NFW3 model at the begin- 
ning of the simulation; the dotted line is the profile at the 
end of the simulation, i.e. after the complete orbital de- 
cay of a model of M54 that is launched in orbit within the 
NFW3 halo. It is interesting to note that the theoretical 
profiles are fairly flat in the considered radial range, in 
good agreement with the data. It may be conceived that 
the (baryonic) nucleus of Sgr lies within the inner cusp 
of the NFW halo of dark matter that is presumed to em- 
bed t he Sgr galaxy (jlbata fc Lewis|[l998l : lMaiewski et al.l 
l2003l and references therein). The potential in this inner 
region would be dominated by the DM cusp that will im- 
pose the flat dispersion curve to the observed Sgr,N stars. 
It does not seem necessary to assume a strict correlation 
between the typical size and/or density proflle of the DM 
cusp/core and those of the embedded stellar nucleus: the 
DM overdensity would have simply provided the "local" 
minimum of the overall potential well to "attract" the 
largest density of the infalling gas that formed Sgr,N. It 
may be interesting to note that the NFW3 model, whose 
velocity dispersion proflle is shown in the lower panel of 
Fig. refsisag, encloses a mass of 3.2 x IO^Mq within the 
range of projected radii covered by our data {vp ~ 70 
pc), and that the region in which its projected density 
proflle is well fitted by a power-law, i.e. the inner cusp, 
has a size of rp < 100 pc, similar to Sgr,N. 

Unfortunately, we have no conclusive observational ev- 
idence in support of the above scenario. Note also that a 
DM halo with a core much larger than the nuclear scale 
would also display a flat dispersion curve, in this region 
(as discussed above for the case of the C=0.90, Vc = 224' 
King model); it is just the small-scale inner cusp as a 
possible seed for the formation of a barionic nucleus that 
can make the NFW option slightly more attractive. 

® We stress here that we are not attaching any particular physical 
meaning to King models while fitting Sgr,N profiles. This is just 
another way to put in evidence the differences between Sgr,N and 
M54 and a convenient reference model to param etrize the observed 
SB profile l|Muno2. Maiewski fc Johnstonll2007l '). 



A flat velocity dispersion proflle may arise as a conse- 
quence of strong tidal di sturbance, in mass-follow-light 
mode ls of dwarf galaxies (jMunoz. Maiewski fc JohnstonI 
I2Q07D . While Sgr is clearly in course of tidal disruption, 
here we are dealing with the kinematics of the innermost 
~ 70 pc of the galaxy, that should (reasonably) be con- 
sidered as untouched by the Galactic tides, since the tidal 
radius of Sgr is as large as a few thousands of parsecs (the 
formal King' s major axis limiting radius is rt > 10 kpc, 
according to lMaiewski et al]|2003D . 

Finally, it may appear as a curious coincidence that 
M54 and Sgr,N have similar spatial scales, if they have in- 
dependent origin. In this sense we can only note that the 
distributions of sizes and luminosities of globular clusters 
and stellar galactic nuclei largely overlap, so the observed 
similarity rnay no t be particularly odd (see Fedcri ci et al.l 
120071: lB6keill2007[ ). Moreover, it should be noted that the 
core and half-light radii of M54 are ~ 2 times larger than 
those of Sgr,N, clearly a non-negligible difference. Fi- 
nally, while the best fltting King models have essentially 
the same C parameter, the two observed proflles appear 
to have significantly different slopes at large radii, where 
the profile of Sgr,N departs from the King model. 

4.3. Rotation 

We searched for signal of rotation in the two samples, 
also taking into acco unt the correction for p erspective 
rotation, according to Ivan de Ven et al.l (pOOo) an d tak- 
ing th e proper motion of Sgr/M54 from Dinescu et al.l 
(|2005( ): in the present case the absolute value of the max- 
imum correction is < 1.0 km/s. Given a rotation of the 
coordinate axes by an angle 9 

7j = Xcos{6) - Ysin{e) 

X = Xsin{9) + Ycos{9) 

we tried all the possible values of 9 in steps of 1° . For each 
adopted 9 we computed the median velocity of stars with 
negative and positive r/. The difference between the two 
median velocities is a robust measure of the amplitude of 
any systemic motion around an axis tilted by 9 degrees 
from the x axis, i.e. of the rotation of stars about this 
axis. 

In some cases, ordered patterns of Vr as a function 
of r] emerged (in particular for the M54 sample), but 
their amplitude was quite weak, typically lower than 2.0 
km/s. While we can clearly exclude the presence of ro- 
tation signals stronger than this, we cannot exclude that 
such weak patterns are in fact real. In Fig. [15] we show 
the best case for the M54 sample compared with the 
end product of one of the N-body simulations that are 
described in Sect. 5, below. We refrain from attaching 
any significance to such a weak pattern, in particular as 
it has an amplitude similar to the typical uncertainty 
of the individual velocity measures and as it is much 
weaker than the typical velocity dispersion of the sam- 
ple. We just want to draw the attention of the reader 
to the fact that the orbital decay of a massive cluster 
within the parent dwarf galaxy seems to naturally pro- 
duce these kinds of weak velocity gradients in the in- 
ner part of t he sinking cluster, at the end of the decay 
process (see iMuiioz. Ma iewski fc JohnstonI l2007l for a 
deeper discussion of the effects of tidal disruption) . How- 
ever, it is quite clear that the detection at r > 10' of stars 
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Fig. 18. — Left Panels: Rotation curve of M54 assuming a rotation axis tilted by 6 = 60° East of North. AVr is the difference between 
the observed radial velocity and the global mean systemic velocity of Sgr,N ((Vr) = km/s). Upper panel: individual stars. Lower panel: 
running mean with bin width of 2' and step of 0(5. Right Panels: Rotation in the remnant of the model of M54 at the end of the simulation 
cll5_b. The coordinate system is chosen such that x,z projects on the plane of the sky while y is along the line of sight and the orbit of the 
cluster is in the x,y plane. Upper panel: l.o.s. velocity as a function of x for the innermost 300 pc. The vertical thin lines enclose a region 
of the same dimension as the one studied here; the large dots are the running mean of the velocity computed with the same bin width and 
step as in the real case. Note the tidal tails emerging from the main body of the remnant a,t x ±50 pc. Lower panel: zoomed view of the 
inner region of the plot. In this case x has been converted in arcmin by placing the remnant at the same distance as the real cluster. The 
diagram has the same scale as the lower left panel for an easy comparison. 



having metallicity compatible with M54, possibly lower 
velocity dispersion than the surrounding metal-rich Sgr 
population, and, above all, having significant differences 
in mean systemic velocity compatible with a tidally in- 
duced rotation like that shown in the upper right panel 
of Fig. [iHl would be the final "smoking-gun" of the or- 
bital decay of M54 to its current position. On the other 
hand the failure to find such a component would not be 
sufficient to rule out the hypothesis as, for instance, the 
density of stars in the tidal tails is expected to be very 
low and the "rotation" signal may be greatly weakened 
by unfavorable orientations of t he orbital plane with re- 
spect to the line of sight (,Munoz. Maiewski fc Johnstonl 
[200l . 

5. THE KINEMATICS OF M54 AND SGR,N: N-BODY 
SIMULATIONS 

M05a used simple analytical formulae to verify the 
plausibility of the hypothesis that the "exact" coinci- 
dence between the positions of M54 and Sgr,N is due to 
the cluster progressively spiraling into the center of den- 
sity of the Sgr galaxy due to dynamical friction. They 
found that the observed status of the globular cluster 
system of Sgr is in full agreement with this hypothesis. 
The orbit of the massive M54 cluster is expected to decay 
completely within one Hubble time if it is born within 
3 core radii {'^ 5 kpc) from the center of Sgr, while the 
other, much less massive Sgr clusters (Ter 7, Ter 8, Arp 2) 
would be essentially unaffected by dynamical friction if 
they were born outside 1 core radius. While promising, 
this first result demands further investigation. One very 
interesting question, for instance, is: does the actual or- 
bital decay of M54 under realistic conditions lead to such 
a nearly perfect coincidence between M54 and Sgr, N as 
observed? To answer questions like this we performed 
the suite of N-body simulations that is described in de- 



tail below. While examining the results of our N-body 
experiments it should be carefully considered that: 

• The main purpose of our simulations is to study the 
effect of dynamical friction in the specific case of 
M54 orbiting within the Sgr galaxy, with particular 
focus on the final status of the system. Detailed 
and more general theoretical analyses of the effects 
of dynamical friction on globul ar clusters within 
a dwa r f galaxy can be fou nd in lOh. Lin fc Richer 
(|2g00f ). lRead"e t al. (2006), Fuiii. Funato fc Making 
(120061). iMiocchi et all (120061). a nd 
[Sinchez-Salchedo. Reves-Iturbide fc Hernanded 
(|2006D . and references therein. 

• We do not intend our simulations to be exhaustive 
of all the possible cases. This would be prohibitive 
since, for example, we do not even know what is 
the total mass of the boun d part of the Sgr galaxy 
fsee iMaiewski et al]|2003l and references therein). 
Our aim is to verify if a broadly realistic case can 
produce the final outcome we actually observe, i.e. 
a cluster that has nearly exactly the same position 
and velocity as the nucleus of its host galaxy. 

• Here we study the case of a massive cluster orbiting 
within a model of the Sgr galaxy which evolves in 
isolation. We verified that this assumption does 
not seriously affect our main results. We repli- 
cated some of our experiments launching the Sgr 
model into a r ealistic (but static) Galac tic potential 
(model 2d of iDehnen fc BinnevI [l998l ): we found 
that if M54 was not lost into the tidal tails of the 
disrupting S gr galaxy (as it happened for some Sgr 
clusters, see iBellazzini. Ferraro fc Ibata 2003aD. it 
plunged to the center of Sgr approximately on the 
same timescales as the isolated models. 
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• A detailed comparison with the existing literature 
on dynamical friction would be clearly beyond 
the scope of the present paper (see, for instance 

IColpi. Mayer fc Governatol Il999l and references 
therein). In the present context, suffice to say 
that our results are in reasonably good agreement 
with the analytic esti mates by M05a as well a s 
with the results by lOh. Lin fc Rich eil (|2000D : 

lSanchez ~Salchcdo . Reyes-Iturbi dc fc Hc rnandej 
(|2006D : lRead et al.l (|2006l . and references therein). 



5.1. General features of the N-body experiments 

The simulations were performed with falcON, a fast 
and momentum-conserying tree-code (Dchnen 200C| 
I20021 ). within the NEMO cnyironment (Tcubcn 1995). 
Grayity was softened with the kernel 'P2' (Dehnen 2000), 
with a softening length of 3 pc (except for two special 
cases, bh_15p and bh_02a, discussed in Sect. 5.2, below). 
The minimum time-step was 1.9 x 10~^ Gyr, and the 
tolerance parameter 9 — 0.6. 

To model the Sgr ga l axy w e adopted a truncated 
iNayarro. Frenk fc Whiti (|1996l hereafter NFW) model 
of the form: 

sech(r/6) 

For the model NFWl, the parameters a = 0.1 kpc, 
h = 2.5 kpc, were selected, which together with a maxi- 
mum circular speed of Vc — 30 km/s, giyes a total mass 
of Mtot = 2.4 X 10® Mq. This rather massiye model 
was chosen t o approximate to the sort of dense system 
suggested by llbata fc Lewis! ()1998| ) that would not be 
rapidly destroyed by Galactic tides. A second, lower 
mass system (NFW2) with a = 0.1 kpc, h = 2.5 kpc 
and Vc = 15 km/s {Mtot — 6.1 x 1O''M0) was also 
simulated. The motiyation for this second model was 
that the central yelocity dispersion in this case should 
be similar to the obseryed dispersion of Sgr (while keep- 
ing a and b identical to the first model) . The NFW3 with 
a = 0.5 kpc, b = 15.0 kpc and Vc = 17 km/s was chosen 
to haye both a size and a central yelocity dispersion sim- 
ilar to the present-day main body of Sgr, irrespectiye of 
its robustness to Galactic tides. 

In all cases the NFW halos were modeled with 10^ 
particles. M54 was modeled as a single massiye particle 
of 1, 1.5 or 2 million solar masses, as if it were a black 
hole (bh simulations, as opposed to cl simulations, see 
Tab. [TOl and Sect. 5.2, below). The particle is launched 
within the NFW halo from (x,y,z)=(2,0,0) kpc or (4,0,0) 
kpc^ (note that the core radius of Sgr is Tc ~ 1.7 kpc 
iMaiewski et al.|[2003D . with yelocities in the y direction 
in the range 4.0 km/s < Vy < 12.0 km/s, while Vx = 
Vy = 0.0 km/s, hence the orbits of the M54 models lie 
in the x,y plane. As a mere conyention to simplify the 
discussion, we consider the y direction coincident with 
the line-of-sight between us and the center of Sgr, and 
the X and z directions as projected on the plane of the 

® In the coordinate system in which (x,y,z)=(0,0,0) kpc corre- 
sponds to the center of density of the model halo and a point having 
(Vx, Vy, Vz) = (0, 0, 0) km/s is at rest with respect to the center of 
the system. The adopted range of initial distances from the center 
of the parent galaxy is typical of the other clusters residing in the 
main body of Sgr (Ter 7, Arp 2, Ter 8) which are located between 
~ 3 kpc and ~ 5 kpc from the center of Sgr and should have orbits 
that are stable against dynamical friction (see M05a and below). 



sky. When needed we will select particles according to 
their projected radius Tp = ^/x'^ + z^. 

The initial conditions of the yarious simulations are re- 
ported in Tab. 1101 The simulations were usually stopped 
when it was clear that the orbit of the M54 particle had 
reached a stable configuration at the center of the host 
halo(ttot)- The decay time {id) is defined as the epoch 
at which the 3-D apocentric distance of the M54 parti- 
cle become rapo < 30 pc. Note that around this limit 
the adopted time step becomes too small with respect to 
the orbital period of the infalling massiye particle, hence 
the eyolution of the orbit cannot be followed further. To 
proyide a quantitatiye idea of the eyolution of the ec- 
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we report the eccentricity at the first (ei) and fifth (55) 
pericentric passages. The eyolution of the orbital radius 
and of the ayerage yelocity of the M54 models is shown 
in Fig. m 

5.2. The evolution of the orbit of M54 models 

From the inspection of Tab.fTOland Fig.[T9]we can draw 
the following conclusions: 

1 . In spite of the wide range of initial conditions that 
has been explored, the orbit of all the M54 models 
decayed completely in less than 3 Gyr (5-15 orbits). 
At the end of our simulations M54 is always yirtu- 
ally at rest (mean yelocity < 1 — 2 km/s) at the 
yery center of the host halo (within a few softening 
lengths). Hence the hypothesis that M54 reached 
its presently obseryed status by dynamical-friction- 
driyen orbital decay appears completely realistic 
and yiable (within the limitations imposed by the 
resolution of our N-body experiments). 

2. Two special simulations (bhl5_p and bh02_a) haye 
been performed at a lower resolution (softening 
length of 10 pc) as the BH was not expected to 
reach the densest regions of the host halo (NFW3). 
bhl5_p explores the case of a M54-like point mass 
decaying from somehow extreme initial conditions 
(large distance - Xq = 5.0 kpc - and high yelocity 
Vy,o = 12.0 km/s). In six Gyr the orbit of the BH 
is sufficiently decayed to bring it within 2 kpc from 
the center of the halo. The other simulations indi- 
cates that from this point it will require less than 3 
Gyr to reach the yery center of the halo. From the 
bhl5_p simulation we can deduce that for starting 
distances larger than 5.0 Kpc the times required 
for a complete decay begin to be comparable to 
the Hubble time, in good agreement with the re- 
sults by M05, and also that it is quite possible that 
M54 haye reached the yery center of Sgr only in 
recent times, if it was born sufficiently far away 
from the center. The analytical computations by 
M05 also indicated that clusters with a mass sim- 
ilar to the other GC residing in the main body of 
Sgr (Ter 7, Arp 2, Ter 8, M < 2. x IO^Mq), would 
haye infinitely long decay times (i.e. stable orbits) 
if they born at r > 1 kpc from the center of Sgr. 
As the closest of these clusters lie at ~ 2.9 kpc 
from the center of Sgr M05a concluded that all of 
them are on stable orbits. Simulation bh02_a is in- 
tended to yerify this conclusion, by letting eyolye a 
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TABLE 10 

N-BODY SIMULATIONS: FUNDAMENTAL PARAMETERS AND RESULTS. 





halo 


Mhalo 






N , 


Xo 


Vy,o 




y^tot 


£1 


^5 






Mq 








kpc km/s 


Gvr 


Gvr 






bh2 a 


NFWl 


2.4 X 10* 


10^ 


2. Ox 10^ 






A n 
4.U 


0.37 


0.84 


0.86 


0.67 


bh2 b 


NFWl 


2.4 X 10* 




2 Ox 10^ 




2.0 


8.0 




1 28 


u. / o 




bh2 c 


NFWl 


2.4 X 10* 


10^ 


2. Ox 10^ 




z.U 


ion 
iz.U 


1.28 


1.58 


0.56 


0.45 


bh2 d 


NFWl 2.4 X 10** 




2. Ox 10*^ 




A n 


A n 


1.22 


1.28 


0.87 


0.65 


bh2 e 


NFWl 


2.4 X 10* 




2. Ox 10^ 




A n 


o.U 


2.90 


3.00 


0.68 


0.65 


bhl_a 


NFWl 


2.4 X 10* 




l.Ox 10*^ 


I 


2.0 


4.0 


0.67 


0.74 


0.83 


0.72 


bhl_b 


NFWl 


2.4 X 10* 




l.Ox lO'' 


I 


2.0 


8.0 


1.19 


1.34 


0.72 


0.58 


bhl_N 


NFW2 


6.1 X 10^ 


10^ 


l.OxlO'' 


1 


2.0 


4.0 


2.25 


2.50 


0.36 


0.13 


cll.2_N 


NFW2 


6.1 X 10^ 


10^ 


1.2x10^* 


10* 


2.0 


4.0 




2.50 






bhl5.a 


NFW3 


4.2 X 10* 


10^ 


1.5x10^ 


1 


2.0 


4.0 


0.80 


1.00 


0.76 


0.71 


bhl5.b 


NFW3 


4.2 X 10* 


IQS 


1.5x10^ 


1 


2.0 


8.0 


1.25 


1.70 


0.61 


0.44 


bhl5_c 


NFW3 


4.2 X 10* 


10^ 


1.5x10*^ 


1 


2.0 


12.0 


2.05 


3.00 


0.40 


0.36 


bhlS.d 


NFW3 4.2 X 10* 


lo--^ 


1.5x10^ 


1 


3.0 


8.0 


2.10 


2.42 


0.60 


0.52 


cll5.b 


NFW3 


4.2 X 10* 


ID--^ 


1.5x10^ 


10* 


2.0 


8.0 




2.00 






bhlO.b 


NFW3 


4.2 X 10* 


lO'"^ 


1.0x10'^ 


1 


2.0 


8.0 


1.60 


2.05 


0.61 


0.49 


bhl5.p 


NFW3 4.2 X 10** 


IQS 


1.5x10^^ 


1 


5.0 


12.0 




6.00 


0.25 


0.15 


bh02_a 


NFW3 


4.2 X 10* 


iqs 


2.0x10* 


1 


2.0 


4.0 




6.00 


0.75 


0.75 



Note. — NFWl; Vc = 30 km/s; NFW2: Vc = 15 km/s. NFWl, NFW2: a = 0.1 and b = 2.5. NFW3: Vc = 17 km/s; a = 0.5 and 
b = 15.0. t^j: time at which the orbit is decayed to a distance from the center r < 30 pc. ttot- total time of the simulation. cll.2_N: t^j 
cannot be defined as in the above cases since the cluster is no more self-bound at t > 0.8 Gyr; at t ~ 1.0 the remnant is approximately at 
the center of the host halo. The case of cll5_b is analogous. 

* At the end of the simulation the BH has reached a distance of ~ 2 kpc from the center of the halo. 

* The orbit of the BH remain stable over the whole duration of the simulations. 
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Fig. 19. — Orbital evolution of the single massive particle representing M54 in the N-body experiments described in Tab. llOl Left panels: 
distance from the center of the host NFW halo as a function of time. Right panels: mean one-component velocity of the massive particle 
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BH as massive as the most massive Sgr cluster (ex- 
cept M54, i.e. Arp 2, M < 2. x ICMq, adopting 
the same {M/L)v of M54, wich seems appropriate 
given the age and metaUicity of the cluster) from 
a starting point Xq = 2.0 kpc and a low initial ve- 
locity [Vvfi = 4.0 km/s). After 6 Gyr no sign of 
decaying is noticed, the orbit is absolutely stable, 
in excellent agreement with the results by M05a. 

The effect of dynamical friction is to progressively 
decrease the radius, the period and the eccentricity 
of the orbit. The mean velocity, on the other hand. 



appears to drop suddenly in the last phases of 
the decay, possibly when the point-mass particle 
reaches the dense central cusp of the host NFW 
halo. This may suggest that the presence of a 
central cusp may be a crucial ingredient to lead 
to the observed phase-space coincidence between 
M 54 and Sgr,N, in particular if the recent results 
bvlSanchez-Salchedo. Reves-Iturbide fc Hernanded 
(|2006D and lRead et all (|2006[) are considered. This 
suggestion is worth being followed up with specific 
high-resolution simulations adopting different 
kinds of host galaxy halo models (such as King 
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Fig. 20. — Evolution of the M54 cluster model within the NFW 
model of the Sgr galaxy, for the simulation cll.2_N (see Tab. llOll . 
In each snapshot of the x,y distribution we plot only the points 
of the M54 model, for clarity. Concentric circles of radius 60 pc 
and 500 pc are plotted in each panel, for reference. In the top 4 
panels we also plotted the orbit of the bhl_N model, to show the 
approximate path of the center of mass of the cluster model. 

models, for instance). Recent studies indicates 
that dynamical friction is probably much less 
efficient in cored than in cusped structures ( see 
Sanchez - Salchedo. Reves-Iturbide fc Hernandea 
2006; Re ad et alJl2006HGoerdt et al.ll2006n . 

To check the impact of our adoption of a point-mass 
model of M54, we repeated the simulation bhl_N and 
bhl5_b adopting a 10^ particle King model resembling 
the real cluster as much as possible {Wq — 8.0, rt = 60 
pc; simulations cll.2_N and cll5_b of Tab. [T0|); in the 
case of cll.2_N the mass was adjusted to M = 1.2 x 
10^ M© to have a central l.o.s. velocity dispersion of ctq ~ 
15 km/s. All the final results are fully consistent with 
those obtained in simulations with a point-mass model 
for M54, for the purposes of the present study, hence we 
do not further distinguish between these two classes of 
models. To illustrate the typical behavior of these live 
M54 models in Fig. [20] we show a series of snapshots of 
the evolution of the cluster model of simulation cll5_b. It 
is interesting to note that in this simulation (as well as in 
the cll.2J>J one) the cluster suffers strong tidal disruption 
from the host halo and the final relic sitting nearly at rest 
at the center of the halo is probably unbound, at odds 
with the real case. The very fact that the actual cluster 
is still bound may provide useful constraints on the mass 
distribution within the Sgr galaxy, that can probably be 
explored with a systematic suite of N-body simulations. 

6. SUMMARY AND DISCUSSION 

We present the results of a large radial velocity sur- 
vey of stars in the nucleus of the Sgr galaxy and in the 
globular cluster M54, which lies at the center of the nu- 
cleus itself. From high signal to noise Keck-DEIMOS 
and VLT-FLAMES spectra, we obtained accurate radial 



velocities (ey^ ~ 2 km/s) and metallicities from the Ca 
triplet for 1152 candidate RGB/RC stars of M54/Sgr,N. 
Selecting by the position in the CMD, by radial velocity 
and by metallicity, we obtained two reliable and clean 
samples of 321 Sgr,N members and 425 M54 stars and 
we used these samples, further cleaned from 3(t velocity 
outliers, to study the kinematics of the two systems. 

In support of this analysis we used archival HST/ACS 
images to refine the estimates by M05a of the structural 
parameters of Sgr,N. 

We also performed a suite of N-body simulations of the 
orbital decay (driven by dynamical friction) of a massive 
cluster into a host dwarf spheroidal galaxy, mimicking 
the M54 - Sgr system. We focused our attention on the 
possibility to realize a phase-space match between the 
infalling cluster and the density cusp of the host NFW 
halo as close as that observed with M54 and Sgr,N. 

The main results of our analyses can be summarized 
as follows: 

1. We have obtained new, much more reliable esti- 
mates of the total luminosity and of the typical 
size of Sgr,N, confirming that it is significantly less 
luminous than M54 and that it has a different lu- 
minosity distribution with respect to the cluster. 

2. The systemic radial velocities of M54 and Sgr,N are 
identical within ^ ±1 km/s. Coupling this result 
with those by M05a, it appears that the two sys- 
tems coincide in phase-space to within the obser- 
vational errors, at least for what concerns position 
in space and radial velocity. 

3. The velocity dispersion profile of M54 is in good 
agreement with the theoretical profile of the King 
model that best fits its surface brightness, at least 
over a range of 30 core radii (r < 3.5'), but possibly 
also beyond this radius (see Sect. 4.2). In particu- 
lar, the velocity dispersion drops from a = 14.2 
km/s at r = 0, to cr ~ 5.3 km/s at r ~ 3.5'. 
The observed velocity profile strongly suggests that 
M54 behaves as an ordinary, self-bound globular 
cluster, at least in its innermost and densest re- 
gion. In our view, this indicates that even if at the 
present epoch the motion of its stars is driven by 
the potential produced by the overall mass distri- 
bution within Sgr,N (see below), its velocity profile 
keeps memory of its original nature, that is an or- 
dinary massive globular cluster orbiting within the 
Sgr galaxy. A turn-over of the profile to cr = 8.8 
km/s around the tidal radius of the cluster may 
suggest that some degree of tidal disruption is oc- 
curring in the outskirts of the cluster; however the 
alternative hypothesis that the sample is contam- 
inated by metal-poor Sgr stars at large radius is 
also viable and seems more likely, at the present 
stage (Sect. 4.2). 

4. The velocity dispersion profile of Sgr,N is consis- 
tent with being flat at cr ~ 10 km/s over the whole 
considered radial range (r < 9'). The profile is 
hardly compatible with any realistic King model 
roughly reproducing the observed surface bright- 
ness profile of Sgr,N. More generally, the fact that 
the velocity profile is flat in a radial range in which 
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the surface density declines by a factor of > 1000 is 
strongl y uncompatible with any mass-follows-light 
model ( Gilmore et al.|[2"007( ). Conversely, a realistic 
NFW model reasonably reproduces the observed 
velocity profile and, in principle, does not conflict 
with the presence of an overdensity of baryons at 
its center, like Sgr,N. However, it has to be recalled 
that the NFW model does not make any definite 
prediction on the SB profile of the embedded stel- 
lar nucleus and the compatibility of the observed 
SB profile of Sgr,N with the NFW model that that 
fits the velocity profile is not established. 

5. M54 and Sgr,N have definitely different kinematical 
properties. In particular, the velocity dispersion 
profiles are very different: in the radial range 1.5' < 
r < 6.5' the statistical significance of the difference 
in the velocity distribution is very large. 

6. We have provided observational evidence that the 
velocity dispersion profile of Sgr remains flat from 
r ~ 1' to r ~ 100' and that there is no apparent 
transition in the velocity proflle corresponding to 
the onset of the stellar nucleus, at r < 10'. This 
fact as well as those listed above strongly suggest 
that Sgr,N and M54 had independent origins, as we 
would expect that, if the cluster provided the mass 
seed to collect the Sgr gas that later formed the 
metal-rich nucleus, Sgr,N stars would have shown 
a declining velocity dispersion proflle, compatible 
with a mass-follows-light distribution. 

7. Our N-body simulations that follow the orbits of 
massive clusters (1-2 xIO^Mq, representing M54) 
within different NFW halos (representing the Sgr 
galaxy) show that for a large range of initial dis- 
tances and relative velocities, the orbit of the clus- 
ter decays completely by dynamical friction within 
3 Gyr, at most. Moreover, at the end of the simu- 
lations, the cluster is perfectly concentric with the 
cusp of the host halo (within the resolution of the 
simulation) and the difference in average velocity 
is always less than ~ 2 km/s. Hence, the observed 
phase-space coincidence between M54 and Sgr,N 
can be naturally explained by the "dynamical fric- 
tion hypothesis" (M05a). 

8. According to FC06 the mass of Central Massive 
Objects, independently- if they are Black Holes or 



stellar nuclei, is 



galaxy (^ 



1000 



of the mass of the host 
0.003). Using our estimates 
for the total mass of Sgr,N and M54 we obtain 
the following estimates for the total mass of the 
Sgr galaxy, depending on whether we assume is 
the M54, Sgr,N or the sum of the two (as we 
would do if we observed the system at the distance 
of Virgo cluster) is the Central Massive Object: 
Msgr = 1.0 X 10^ Mq, Msgr = 2.1 X 10^ Mq, and 
Msgr = 3.0 X 10^ Mq, respectively. These values 
are in very g ood agreement w i th the independent 
estimates by iMaiewski et al.l ()2003[ ) that ranges 
from Msgr 5.8x 10* Mq to Msgr = 6.9x 10^ Mq. 
Hence both Sgr,N and M54 singularly or taken to- 
gether have a mass compatible with being the CMO 
of the Sgr galaxy. 



9. Both Sgr,N and M54, as well as the combination 
of the two, when placed in the My vs. logr/j di- 
agram lie in a region that is populated by globu- 
lar clusters and galactic stellar nuclei. They are 
also compatible with the Color-Magnitude relation 
of nuclei shown by FC06 and with the My vs. a 
relation sat isfled by globular clusters, nuclei and 
UCDs (see iGeha. Guhath akurta &: van der Marell 
12001 lEvstigneeva et al. 200^ as shown in Fig. \TT\ 
Hence the structure and dynamics of Sgr,N, M54 
and their combinations are fully compatible with 
other galactic nuclei. 

10. A detailed study of the mass profile of Sgr,N 
and M54 using the Schwarzschild ni ethod (see 
iRix et al.l [19971 : Ivan de Ven et al.ll2006l . and refer- 
ences therein) is currently ongoing and the results 
will be presented elsewhere (Ibata et al., in prepa- 
ration) . 

These findings lend very strong support to the scenario 
proposed by M05a to explain the M54/Sgr,N system: the 
nucleus of the galaxy formed in situ, at the bottom of 
the potential well of the Sagittarius galaxy; the globular 
cluster M54 was independently driven to the same site 
by dynamical friction. As a complement of the above 
conclusions, it must be recalled that the stellar popula- 
tion that dominates Sgr,N formed s everal Gyrs later than 
M54 (see [Montcgriffo et al. 199 8t ISaraiedini fc LavdenI 
1997; Lavden fc Saraiediniii2000l : lBeriazzini et al.ll2006al : 
Siegcl et al. 2003, and references therein) . In the present 
context this does not appear particularly relevant, but it 
should be kept in mind that M54 can have reached the 
very center of Sgr after Sgr,N forrned its s tars (within the 
last ~ 5 - 9 Gvr IBellazz ini et al.l l20"06al) . or even before 
this, depending on the birthplace of the cluster within 
the Sgr galaxy. In any case it is very likely that both 
the cluster and the processed Sgr gas were independently 
driven - by different mechanisms - to the bottom of the 
potential well of the Sgr galaxy, i.e. to the center of 
its Dark Matter halo. A conclusive proof that M54 was 
driven to its present position by dynamical friction could 
be provided by the successful detection of genuine extra- 
tidal stars at large distances from the cluster center as 
envisaged in Sect. 4.3, while, as said their non-detection 
will not disprove the above scenario. 

As a possible alternative to this view (or as an ex- 
treme version of it) it can be conceived that M54 formed 
at the bottom of the overall potential well of the Sgr 
galaxy since the beginning. In this framework M54 and 
Sgr,N are just the results of two subsequent episodes 
of star formation both occurring at the very center of 
Sgr, the second from enriched gas that was infalling 
on very different orbits with respect to the gas that 
formed M54, thus resulting in different final stellar kine- 
matics. We regard this possibility as much less likely 
with respect to the "dynamical-friction" scenario de- 
picted above. First, if M54 is considered as part of the 
(metal-poor) field population of Sgr, its presence would 
be at odd with the global metal licity gradient o bserved 
by .A lard (2001), BeUazzini. Fcrra ro fc Buonannol (fl999). 
ISarajedini fc LaydenI (|1997l) and others in Sgr, as well in 
all dwarf spheroidal galaxy studied to date (more metal 
rich and younger populations are preferentially found at 
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Fig. 21 . — Ssr.N and M54 (fill ed triangles) in the Mv vs. a 
plane (see 'Evstigneeva et al."2003). Both systems lie on the locus 
common to globular clusters, d E nuclei, Dwarf-Globular Transi- 
tion Objects (DGTO Hasegan c t al.l 120051) , and, possibly, UCDs. 
The same is true for their combination (Sgr,N+M54, open trian- 
gle), i.e. for a system having luminosity equal to the sum of the 
luminosities of M54 and Sgr,N and having a as estimated from an 
integrated spectrum obtained with a 1" slit from the distance of 
the Virgo cluster of g alaxies (ct ~ 12 A km/s). Data for Galac- 
tic globu lars are from iMackev fc van den Berghl | |2005| , My) and 
IPrvor fc Mcylan ( 1993, ct); for M31 globulars we^took integrated 
magnitudes from the R evised Bologna Catalog ( Galleti et al. 2004) 
and CT fr om Diorgovski ct al. (1997); data for dE and dE,N arc from 
IGeha, G uhathakurta fc van dcr Marel (2002, 2003); UCDs data 
are from Drinkwater et al.| 820031) and Evstigneeva ct al. (2007); 
data for DGTO are from lHasegan et a l. (2005). Some remarkable 
brigh t clusters have been labeled, for reference (see iFederici et al.l 
[20071) . 

the center of dSph's lHarbeck et al.|[200ll ). While dE nu- 
clei have been indicated as a possible exceptions to this 
trend, the ingestion of large metal poor globular clus- 
ters in pre-existing metal rich nuclei seems one possible 
natural way to reconcile the generally observed metal- 
licity gradients and the presence of nuclei that are blue r 
than their parent galaxy ()Lotz. Miller fc Fe rguson"2004^ . 
Second, as both M54 and Sgr,N would have formed from 
gas falling into the same potential well, the reason for the 
different kinematics remains to be explaned, while it is a 
natural outcome if M54 formed elsewhere as a classical 
globular cluster. 

6.1. The process of galaxy nucleation 

If we take for demonstrated the above concluding re- 
marks, we can ask what we have learned about the pro- 
cess of galaxy nucleation from the case studied here. 
Concerning the two main mechanisms that have been 
put forward in the literature, i.e. (a) formation of the 
nucleus by infall of globular cluster (s) to the center of 
the galaxy, or (b) in situ formation by accumulation of 
gas at the center of the potential well and its subsequent 
conversion into a stellar overdensity (see Sect. 1. and 
iGrant. Kuipers. fc Phillip^ l2005Hc'6te et al.l[200l . the 
main conclusion that can be drawn from the case of Sgr 
is that both channels are viable and actually both have 
been at work "simultaneously" in Sgr. 

The present analysis has shown that a stellar nucleus 
made of the typical material that dominates the baryonic 



mass budget of Sgr is present in this galaxy, independent 
of M54, as it display the same flat velocity dispersion 
profile as the whole core of Sgr, much different from that 
of the cluster. 

In a likely scenario, the enriched gas from previous 
generations of Sgr stars accumulated at the bottom of 
the overall potential well of the galaxy, until star forma- 
tion transformed it into a stellar nucleus whose surface 
brightness is > 100 times larger than in the surround- 
ing Sgr core, a substructure within a larger galaxy. On 
the other hand, M54 is a (relatively) ordinary massive, 
old and metal poor globular cluster. Independent of its 
birthplace within the early Sgr galaxy, its mass and the 
density of the surrounding medium of the host galaxy 
drove it to the bottom of the Sgr potential well, by dy- 
namical friction. During its trip to the densest central 
region of Sgr, the dense cluster managed to survive the 
tidal force of the host galaxy, hence it reached the present 
position as a (partially?) self-bound stellar system. 

While the mass budget at the center of Sgr is probably 
dominated by the underlying DM halo, M54 dominates 
the overall light distribution : an observer taking pho- 
tometry and/or spectra of the unresolved nucleus of Sgr 
from a distant galaxy (say, a galaxy in the Virgo clus- 
ter) would find that the object looks like a bright and 
blue globular cluster; the integrated velocity dispersion 
would not show any peculiar feature revealing the com- 
posite nature of the observed nucleus (see Fig. [211 above). 
Probably, it would be impossible to disentangle the two 
systems from the integrated light 

Finally, the Sgr case seems to support the observed 
ubiquity of nuclei (see Sect. 1. and C06). The Sgr 
galaxy was able to form a sizable nucleus "twice" and 
with two different formation channels: if either of the 
two channels had not been viable for some reason, the 
galaxy would have ended up with a nucleus in any case. 
The fact that Sgr is the only case of galaxy with a clear 
nucleus among those classified as dwarf spheroidals in 
the Local Group may suggest that the progenitor of 
Sgr was in fact a brighter dE or disc galaxy that has 
been transform ed into a dSph by the interaction with 
the Milky Way (|Maiewski et al.ll2003HMaver et al.ll2007l 
M05b, and references therein) 

6.2. Suggestions for further investigations 

The results presented in this paper suggest several in- 
teresting lines of research that we did not follow up for 
practical reasons. However we feel that it is worth briefly 
mentioning some of them here, as a possible starting 
point for future studies. 

• The result s presented b v iRead et al.l 
(12001. iGoerdt et al.l (120061) a nd 

[Sdnchez-Salchedo. Reves-Iturbide fc Hernandez! 
(120061) . suggest that the complete decay of a 
massive cluster to the very center of the host 
galaxy is much easier and faster if a central 
density cusp is present. It is even possible that 
a central cusp is actually required to bring a 

It would be interesting to check if there is some spectral fea- 
ture that may reveal the composite nature of the "system" , in the 
present case. We plan to do this in the future by combining prop- 
erly scaled synthetic spectra representing the light output of M54 
and Sgr,N. 
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cluster down to the very bottom of the overall 
potential well. If this is the case, the position of 
M54 within Sgr,N would support the existence 
of a central cusp in actual DM halos, a point 
that has been qu estioned by several authors ( see 
Sanchez- Salchedo. Reves-Iturbide fc Hernandez! 



20061 and references therein). The "complete 
infall" of a massive cluster may need a NFW 
cusp and, simultaneously, it may transform the 
cusp into a core by transferring orbital energy 
and momentum to the surrounding "medium" , 
thus possibly providing a self-regulating mech- 
anism that simultaneously prevent the fur- 
ther decay of other clusters (that, in general 
cases, would be quite di fficult, given the ex- 
pecte d decay tim es, see iHernande z fc Gilmord 
19981: lOh. Lin fc Richer 2000; Read c t al.l l2006t 
Sanchez- Salchedo. Reves-Iturbide fc Hernandez! 



20061: iGoerdt et al .I I2006I . and references therein). 
In this conte xt, it may be worth to recall also 
the work by IStrigari et al.l (|2006( ). whose results 
militates against the prese nce of a large-size core 
i n the Fornax dSph, and bv lBovlan-Kolchin fc Mai 
()2004f ) on the resilience of cuspy halos in major 
mergers. These ideas seems worthy of detailed 
theoretical follow up. 

• It has been suggested several times that bright 
and metal poor globul ar clusters may be of 
cosmological origi i i (see iBrodie fc Straded l2006t 
iMaschenko fc Silla l2005l and references therein) . 
If this is the case they may be very intimately 
linked to the earliest phases of formation of galax- 
ies and there may have been plenty of oppor- 
tunities for most of them to become the nu- 
clei of some dwarf galaxy. The tidal stress that 
they suffered during their infall to the center 
of their host galaxies may be at the origin of 
the larger half-light sizes that are observed in 
those of them that have been suggested as pos- 
sible nuclear remnants of ancient dwarf galaxies 
([Freeman fc Bland-Hawthorn' '2002': jFtederici et al.l 
|2007; Mackey fc van den Bcrgh 200l). This kind 
of scenario may be explored with dedicated N-body 
simulations, possibly including gas and stars. It 
would be interesting also to consider in detail the 
results presented here in relation with the scenario 
for the origin of globular clusters recently discussed 
bv lB5kei!(j200'7l) . 

• All the analysis presented in this paper have been 
performed within the standard Newtonian gravita- 
tion theory and dynamics. It may be worth to con- 
sider the observational scenario emerged from this 
study also in the framework of Modi fied Newtonian 
Dynamics paradigms (MOND, see Mil groml l2008l : 
iSanders fc McGaughll2002i . and references therein) , 
even if it may not necessarily be an ideal case. The 
transition between ordinary Newtonian regime and 
MOND regime occurs around r ~ 4.5' — 6.5', de- 
pending on the actual stellar mass of M54-|-Sgr,N, 
that is in a range covered by our data. 

• There is little doubt that the final fate of the Sgr 
galaxy will be its complete tidal disruption. Once 



the large scale stellar body of Sgr will be com- 
pletely dispersed into the Galactic halo, the final 
remnant of this (once) relatively large galaxy would 
be a faint nucleus embedding a bright globular clus- 
ter. An observer lacking any knowledge of the ori- 
gin of this object would conclude that it is a very 
bright and peculiar globular cluster, dominated by 
a metal poor population {[Fe/H] ^ —1.5, possi- 
bly with some spread) but also including a small 
fraction 10%) of metal rich stars (with aver- 
age [Fe/H] ^ —0.4). Also the abundance patterns 
would appear different: for instance, the metal 
poor (M54 ) stars would appear as modera tely a- 
enhanced (|Brown. Wallerstein fc Gonzaled [r999), 
while metal rich stars (Sgr,N ) would h ave solar 
or sub-solar [a/Fe] ratios (Monaco et al. 2005b). 
The radial velocities would reveal that the metal 
rich and metal poor stars have the same sys- 
temic velocity, but different velocity dispersion pro- 
files and slightly different density profiles, possi- 
bly with some rotation in the metal poor com- 
ponent. However the dominance of metal-poor 
(M54) stars would be so high that there would 
be no hint of the presence of a Dark Matter 
component. The half-light radius of the system 
would appear sli ghtly larger with respect to or- 
dinary globulars (iMackev fc van den Berghl [20051 : 
iFederici et al.l l2007l ). Most of these likelv charac- 
teristics of the future remnant of Sgr seem to have 
a counterpart in the widely studied and mysterious 
stellar system m Centauri (se e Norris et al. 199j; 
van Lceuwcn, Hug es fc Piottci|[2 002: Pan cino et al] 
2000- ,2002 . ,2003i : 'lvan de Yen et al. 20061. and ref- 



erences therein), that has been proposed as the 
possible remnant of a nucleated dwarf elliptical 
since Freeman (1993). While there are also notice- 
able differences between the two cases, the analogy 
seems very intriguing and potentially powerful^^. 
It is possible that at least some of the observational 
features of u Centauri that appear so difficult to ex- 
plain may find their natural place within a scenario 
like the one described above. 
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It is interesting to recall tliat w Cen sliares with Sgr some re- 
markable chemical peculiarities. In particular, stars of comparable 
metallicity in the two systems are strongly enhanced in s-process 
elements and the ratio of heavy s-process to light s-process ele- 
ments [hs/ls] is very similar. Furthermore, u> Cen and Sgr are, up 
to now, the only stellar systems known to have deficient [Cu/Fe] 
ratios with respect to the trend i n the Galactic disk and halo 
IIMcWilliam fc Smecker-Hanell2005h 
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